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Abstract 

We review the quantum version of the linear Boltzmann equation, which describes in a non-perturbative fashion, by 
means of scattering theory, how the quantum motion of a single test particle is affected by collisions with an ideal 
background gas. A heuristic derivation of this Lindblad master equation is presented, based on the requirement of 
translation-covariance and on the relation to the classical linear Boltzmann equation. After analyzing its general 
symmetry properties and the associated relaxation dynamics, we discuss a quantum Monte Carlo method for its 
numerical solution. We then review important limiting forms of the quantum linear Boltzmann equation, such as 
the case of quantum Brownian motion and pure coUisional decoherence, as well as the application to matter wave 
optics. Finally, we point to the incorporation of quantum degeneracies and self-interactions in the gas by relating 
the equation to the dynamic structure factor of the ambient medium, and we provide an extension of the equation 
to include internal degrees of freedom. 
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1. Introduction 

The Boltzmann equation, as the basic equation of kinetic gas theory, belongs to the foundations of clas- 
sical statistical mechanics. Apart from its great relevance for concrete applications, it is a paradigm for 
the description of classical systems far from the thermodynamic equilibrium and for the understanding of 
irreversibility and equilibration [1|2|3|4|5] . A closely related equation is used if one wants to describe how 
the motion of a single distinguished test particle is affected by elastic collisions with an ideal, stationary 
background gas. In both situations one deals with the evolution of a probability distribution defined on 
the phase space of a single particle. In case of the Boltzmann equation this is the marginal one-particle 
distribution of a dilute gas, whose self-interaction gives rise to the nonlinear form of the integro-differential 
equation. There is no such self-interaction if an individual test particle is considered, such that the corre- 
sponding equation is linear, like the Liouville equation. Apart form this, the linear equation shares much of 
its structure and many basic properties with the original Boltzmann equation, and it is derived the same way 
using the Stosszahlansatz — it is thus aptly called the linear Boltzmann equation [5]. This equation for the 
motion of a distinguished test particle has found important applications, especially in the study of transport 
phenomena [617] . 

The present review is devoted to the quantum analog of the linear Boltzmann equation. It is therefore 
important not to confuse its classical counterpart with a "linearized Boltzmann equation" . The latter is an 
approximation of the original Boltzmann equation for a self-interacting gas, which is valid when the gas 
is close to equilibrium and the interactions are a weak perturbation. In the linear Boltzmann equation, in 
contrast, the coupling to the gas can be strong and the distinguished test particle may be very far from an 
equilibrium state. The background gas, on the other hand, is taken to be in a stationary state, typically at 
thermal equilibrium, which remains undisturbed by the presence of the test particle. 

In the quantum physics literature a large number of master equations has been termed "quantum Boltz- 
mann equation" . The name is likely to appear whenever the dynamics is amenable to a kinetic description 
based on interaction events that can be described as scattering processes. A first basic result on a quantum 
extension of the original Boltzmann equation was the introduction of a correction accounting for the effects 
of quantum statistics in the classical collision term |8l9ll0j . Moreover, Boltzmann-like kinetic equations ap- 
pear in a wealth of situations with the aim of describing relaxation and transport properties in a quantum 
setting. Most recently, the study of quantum kinetic equations required for the description of degenerate 
Bose gases has been the object of extensive research, see, e.g., the series of papers |llll2ll3ll4ll5ll6fT7] and 
references therein. 

In the field of mathematical physics, major efforts have been undertaken in order to rigorously derive a 
quantum version of the nonlinear Boltzmann equation. Here, one typically starts out from the many-particle 
Hamiltonian with a simple interaction potential and seeks to obtain a nonlinear equation of motion for 
the reduced single-particle state by considering approximations for the hierarchy of equations of motion of 
many-particle distribution functions, see |18I19I20I21I22I23] for reviews on recent results. Similar attempts 
have been devoted to the derivation of the quantum linear Boltzmann equation, using a random potential 
approach, see |24|25|26j for recent results. However, it should be emphasized that even the classical Boltz- 
mann equation has not yet been properly derived in the rigorous sense of mathematical physics. The problem 
is substantially more difficult in the quantum formulation, due to the need for accommodating the unitarity 
of the underlying many-body dynamics and due to the quantum restrictions of the associated phase-space 
representation. 

The situation is more tractable if one asks for the proper quantum version of the linear Boltzmann 
equation. Due to the inherent linearity of the problem, the quantum equation of motion must be given 
by a linear mapping of the statistical operator representing the motional state of the test particle. It can 
be expressed in the framework of quantum dynamical semigroups [27128129] since the use of a Boltzmann 
equation implies that the time evolution is Markovian. It follows that the generator of the irreversible time 
evolution must have the Gorini-Kossakowski-Sudarshan-Lindblad form [30|31j : this guarantees complete 
positivity, thus ensuring that the motion of the reduced single particle state is compatible with an underlying 
many particle state that evolves in a unitary fashion. 
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Fig. 1.1. (Color online) The quantum linear Boltzmann equation describes the motion of a single, distinguished test particle 
in an ambient medium, taken to be an ideal, stationary, and uniform background gas. In the framework of the Markov 
approximation, it is valid for arbitrary initial states of motion, such as the non-classical spatial superpositions found in an 
interferometer. As such it describes the interplay of decoherence phenomena at short time scales with the dissipative effects 
and the thermalizing behavior expected in the longer-term. The interaction with the gas is accounted for in a non-perturbative 
and microscopic fashion, by means of the exact two-body scattering amplitudes. The equation ceases to be applicable if the 
neglect of three-body collisions or correlated scattering events cannot be justified, as in the case of liquids. 

Using the framework of Lindblad master equations, the natural quantum counterpart of the classical linear 
Boltzmann equation can indeed be obtained as a unique master equation. It thus describes how the free 
quantum motion of a test particle, possibly in a very non-classical, delocalizcd state, is modified by the 
presence of a background gas. see Fig. 11.11 As with the classical linear Boltzmann equation, the quantum 
master equation applies provided (i) the background gas is sufficiently dilute, such that three-body collisions 
can be neglected, (ii) the particle interaction is sufficiently short-ranged, such that scattering theory can 
be applied, and (iii) the gas temperature is sufficiently high, such that effects of quantum statistics can 
be disregarded and the Markov assumption is justified. Both the classical and the quantum formulation 
incorporate the two-body interaction between test particle and a gas molecule in a non-perturbative fashion 
by means of scattering theory. While the classical equation contains the differential cross-section to account 
for the molecular collisions, the quantum master equation involves pairs of scattering amplitudes. These 
objects, which feature a complex phase, show up as operator-valued quantities if the master equation is 
stated in a representation independent form. 

The main significance of the quantum linear Boltzmann equation lies in the fact that it serves to describe 
the interplay of two distinct gas-induced phenomena, typically occurring on different time scales and usually 
requiring a non-perturbative description of the coupling to the gas. On one hand, this is the coUisional 
decoherence effect, taking place at short times when the energy exchange between test particle and medium 
can be disregarded. On the other hand, the test particle experiences friction, leading on a longer time scale to 
a slowing down of the particle and its eventual thermalization. The equation should thus predict how a test 
particle, starting out in a delocalizcd state, gradually loses its ability to display quantum interference, due 
to the dissemination of position information effected by the gas collisions. One then expects the quantum 
particle, whose motional state is increasingly hard to distinguish from a classical one, to display the relaxation 
and equilibration behavior predicted by the classical linear Boltzmann equation. It is thus an important 
consistency requirement that the quantum linear Boltzmann equation indeed reduces to the classical version 
once the quantum state of motion is indistinguishable from a classical state. Also in the other relevant 
limiting cases the equation turns into well-established forms, in particular in the limit of quantum Brownian 
motion, of pure coUisional decoherence, and of coherent forward scattering, as will be discussed below. 

In the present review we do not aim at providing a microscopic derivation of the quantum linear Boltzmann 
equation; such a calculation was given recently in some detail in [32j . Instead, we will motivate the form of the 
equation in a more succinct line of reasoning by combining natural physical and mathematical requirements. 
It builds on the structure of the classical equation and on Holevo's characterization of general translation- 
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covariant quantum dynamical semigroups |33| . by exploiting both the required limiting expression of the 
equation and the translational symmetry expected in a homogeneous medium. Apart from the operator 
form, we will also discuss the equation and its limiting cases in different representations, each admitting 
specific physical insights. 

The physics described by the quantum linear Boltzmann equation is discussed and applied in quite dif- 
ferent fields of physics, ranging from mathematical physics and statistical mechanics to matter wave optics. 
Traditionally, quantum versions of the linear Boltzmann equation were introduced in order to provide a 
microscopic explanation of transport phenomena, or to obtain a derivation of quantum Brownian motion. 
These are situations where the relevant quantum states are close to the state of motion of a classical particle. 

Only recently it has become possible to study the influence of a gas environment on the motion of a 
molecule which is very far from a classical state, such as in a superposition of macroscopically distinct 
spatial positions [34j . This is a result of the experimental advances in the interferometry of massive particles 
[35l36j . allowing one to observe the phenomenon of decoherence independently from dissipative effects. 
In fact, the main experimental applications and confirmations of the quantum linear Boltzmann equation 
are found in situations where the quantum coherences described by the statistical operator play a crucial 
role. This applies to the study of the refractive index for matter waves traveling through a gas |37|38j . 
and in particular to the analysis of coUisional decoherence for massive particles |39j . It is worth noting 
that the seminal paper on decoherence by Joos and Zeh }40| . seeking to explain the absence of quantum 
delocalization in a dust particle by the scattering of photons and air molecules, derived and studied what the 
authors called a Boltzmann-type master equation. Two decades later, the long quest for the characterization 
of the phenomenon of coUisional decoherence has now reached a mature theoretical description, permitting 
its quantitative experimental confirmation. 

The article is organized as follows. In Sect. [5] we introduce the equation both in its operator form and 
in the momentum representation. Section [3] provides an intuitive line of derivation, based on reconciling 
the expression of the classical linear Boltzmann equation with the general form of a translation-covariant 
generator of a quantum dynamical semigroup. Various properties and features of the quantum linear Boltz- 
mann equation are then discussed in the following sections, by first addressing phenomena with a classical 
counterpart, and then moving to problems which can only be considered in a quantum setting. Specifically, 
Sect. |4] considers the symmetry properties of the equation and its relaxation and equilibration behavior, 
while Sect. [5] is devoted to the limiting case of quantum Brownian motion, obtained if the test particle is 
much more massive than the gas particles. Section |6] deals with a distinguished quantum effect, the index 
of refraction experienced by a matter wave due to the presence of a background gas, which depends on the 
coherent part of the quantum linear Boltzmann equation. The phenomenon of decoherence is then addressed 
in Sect. [71 both regarding coUisional decoherence and more general forms. Section [S] discusses an alternative 
expression of the quantum linear Boltzmann equation in terms of the dynamic structure factor of the gas, 
suggesting a generalization to the case of quantum degenerate background gases. It further provides the ex- 
tension of the quantum linear Boltzmann equation to include internal degrees of freedom. We finally present 
our conclusions in Sect. [9l where we also point to possible future developments. 



2. The quantum linear Boltzmann equation 

We start by introducing the general expression of the quantum linear Boltzmann equation in its representation- 
independent operator form. In order to understand why the scattering amplitudes and distribution functions 
in this expression must be operator- valued, and why the momentum projections parallel and perpendicular 
to the momentum exchange play different roles, it will be helpful to represent it in the appropriate coordi- 
nates. We therefore proceed by presenting the equation in the momentum basis, which will also facilitate 
the comparison with the classical linear Boltzmann equation. We conclude this section with a discussion of 
the Born approximation and a brief outline of existing approaches for deriving the equation. 
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2.1. The operator form 



The quantum linear Boltzmann equation is a Markovian master equation for the statistical operator p 
describing the motion of a distinguished test particle in a gas. It has the form 

= r^[Ho + ffn(P),p]+/:p, (2.1) 

where Hq — P^/ (2A/) is the kinetic energy of the test particle and (P) describes an energy shift due 
to the presence of the background gas, as defined in Eq. (|2.5p . Here and in the following we use sans-serif 
symbols to denote operators in a single particle Hilbcrt space; variables referring to the test particle will be 
set in upper-case letters. 

The superoperator C in (|2.ip is a linear mapping which accounts for the incoherent part of the coUisional 
interaction with the gas. It can be expressed in Lindblad form as |41| 

Cp = J dQ dk^ [c^Q-V^^L (ki, P; Q) pit (k^^ p. Q) e-^Q X/?. 

(2.2) 



i{p,Lt (ki,P;Q)L(ki,P;Q)} 



where X and P are the position and the momentum operator of the test particle, and the curly brackets denote 
the anti-commutator. The integration is over all momentum transfers Q, and over gas particle momenta p 
from the plane perpendicular to Q, i.e., from = |p G M.^ : p • Q = O}. The function L, which appears 
operator- valued in (|2.2p . contains all the details of the collisional interaction with the gas. It is defined by 



L (p, P; Q) = ( rcl (p^g, P^q) - |, rel (pj_q, P^q) + § 




X\/aMp^Q + — ? + mP|Iq)' (2-3) 



where m and M are the masses of the gas particles and of the test particle, respectively, and = 
mM / (m + M) denotes the reduced mass. It involves also the gas number density rigas, the distribution func- 
tion p, (p) of the gas momenta, and, most importantly, the elastic scattering amplitude / (py, p^) determined 
by the two-body interaction between gas and test particle. The subscripts ||Q and _L Q denote the compo- 
nent of a vector parallel and perpendicular to the momentum transfer Q, such that P||q = (P • Q) Q/Q^ 
and P±Q = P — P||Q respectively. Finally, relative momenta are stated by means of the abbreviation 

rel(p,P)^^p-!i;^P. (2.4) 

TO M 

The Hamiltonian correction in (|2.ip is defined by the function [35] 

H,,{P) = -27rh''^ /dpM(p)Rc[/(rel(p,P),rel(p,P))]. (2.5) 
J 

As discussed in Sect. [51 this energy shift is related to forward scattering. 

The quantum linear Boltzmann equation is valid for non-degenerate, uniform gases with arbitrary sta- 
tionary momentum distribution p (p). However, for concrete calculations we will often tacitly take p (p) to 
be the standard Maxwell-Boltzmann expression 

M/3(P)= 3/2 3 ( I ' (2-6) 



with pp = ■\/2m//3 the most probable momentum at temperature T ~ \/ (k-Q/S). 
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2.2. Momentum representation 

To gain some physical insight into the equation it is convenient to express it in the basis of improper 
momentum eigenstates, P|P) = P|P), which will also help in establishing the relation to its classical coun- 
terpart. In the momentum representation Eq. (j2.2|) can be brought into the form 



(PirpIP') = A/i„ (P, P'; Q) (P - QIpIP' - Q> - ^ (P) + ^^^it (?')] (P|p|P'>- (2-7) 

The "gain term", i.e., the first term on the right-hand side, is here given by 

Mn (P, P'; Q) = / dk^ L (kx, P - Q; Q) L* (ki, P' - Q; Q) (2.8) 



with L defined in Eq. (|2.3p . Importantly, (|2.8p reduces to the rate known from the classical linear Boltzmann 
equation if one approaches the diagonals P = P', 

Min(P,P;Q)=A/rn'(P;Q). (2.9) 
This classical gain rate, which involves the quantum mechanically defined differential cross-section cr (p j , ) = 



l/(P/,P^)l^ reads as 



m Q m . 



(P + Q; Q) = ^ j^^ dk, . + + -P„c, 

xa(^rcl(kx,PxQ)-^,rel(kx,P±Q) + ^) . (2.10) 
It has a natural meaning in the sense that the transition rate 

M^' (P. ^ Vf) ^ (Pf, Pf - P,) (2.11) 

provides the rate of collisions which change the test particle momentum from P^ to P^, or equivalently, 
which lead to a final momentum P f due to a gain of momentum Q = Py — P^ . 

The loss term in (j2.7p is given by the arithmetic mean of the loss rates which appear in the corresponding 
classical equation. The latter are related to the classical gain rate by 

M^i, (P) = 1 dQ Mi (P + Q; Q) . (2.12) 

This relationship, which ensures that the probability is conserved in the classical case, follows directly from 
the structure of (j2.2[) . and it guarantees the normalization of the statistical operator p. 
If fact, by considering Eq. (|2.7p only on the diagonal, 

(P|£p|P) = J dQMi (P; Q) (P - QIpIP - Q) - (P) (P|p|P), (2.13) 

one immediately recovers one of the possible expressions of the classical linear Boltzmann equation, as we 
shall see in Sect. |31 It takes the explicit form 

(P|£p|P) = ^ J ^ |^dk^a(^rel(k^,P^Q)-|,rel(k^,P^Q) + |) 



(2.14) 
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2.3. Born approximation 



The weak coupling limit of the quantum linear Boltzmann equation is obtained if one replaces the scatter- 
ing amplitudes in Eq. (j2.3[) by their Born approximation. This simplifies the equation considerably, since the 
Born amplitude associated to the interaction potential V (x) depends only on the difference of the momenta 
according to 



fs {pf,Pi) = -'i7T^hm^{pf\V (x) \pi 



dx V (x) exp 



27rft2 
^fsiPf-Pi). 

This removes the operator- valuedness of the scattering amplitudes in 
to the form 



because the function 



i(p,P;Q) 



■Ib (-Q) Wm p±q 



m Q 

m, 2 



M 



HQ ' 



(2.15) 
reduces 

(2.16) 



such that the dk^-integration in Eq. (|2.2p can be carried out if one assumes the Maxwell-Boltzmann distri- 
bution (j2.6p . The resulting equation then reduces to the one proposed in [42| . Introducing the operators 



Lb iP;Q) = 



I Pm Hgi^aB (Q) 



exp 



with as = l/sl^j it takes the simpler form 

Cbp = JdQ e*«></''LB(P;Q)p4(P;Q)e-*'^''^'"'-^{p,4(P;Q)is(P;Q)} 



(2.17) 



(2.18) 



where only the integration over the momentum transfers is left. 

We note that Eq. (|2.2p can be simplified also by other limiting procedures. The minimally extended 
Caldeira-Leggett master equation is obtained in the diffusive limit described in Sect. 15. 2[ while the master 
equation of pure coUisional decoherence discussed in Sect. I7.11 is gained in the limit of a very massive tracer 
particle, m/M 0. 



2.4. Survey and comparison of derivations in the literature 



The quantum linear Boltzmann equation has a long history, even though the general form Eq. (|2.2p was 
obtained only recently. The key difficulty lies in the necessity to account for the coUisional dynamics in 
a microscopically realistic and non-perturbative way, leading to a dynamics which is consistent with the 
semigroup structure of a proper Markovian description and which displays the thermalization behavior 
known from the classical description. 

A Boltzmann-like master equation already appeared in the seminal work by Joos and Zeh [40] which used 
a scattering theory description to treat the decoherence effect of environmental collisions on the motion 
of a quantum particle. Despite its similar name the master equation presented there is quite different from 
Eq. (j2.2p in that it linearizes the interaction effect and neglects dissipation, thus leading to an infinite growth 
of the kinetic energy of the test particle on long time scales. Indeed, the equation is meant to describe only the 
short-time loss of coherence in the motional state. However, due to the linearization involved the predicted 
localization rate grows above all bounds for increasingly delocalized quantum states, which may quickly lead 
to unrealistically large decoherence rates. 

This result was improved by Gallis and Fleming [43], who proposed a master equation using the full 
scattering matrix. This equation predicts a saturation of the localization rate at large distances, as one 
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expects for local interactions with the environment, though at a numerically incorrect value. Like the master 
equation by Joos and Zeh, this result does not account for dissipation phenomena. 

To our knowledge, the first proposal pointing to a quantum linear Boltzmann equation which also describes 
energy relaxation and therefore grants the existence of a stationary solution, was given by Diosi [44]. As 
discussed in Sect. 13.4] this result is not equivalent to what we consider the correct quantum linear Boltzmann 
equation, although the diagonal matrix elements in the momentum representation coincide, yielding the 
expected classical equation. 

The correct result including both decoherence and dissipative effects, though restricted to a perturbative 
scattering cross-section evaluated in Born approximation, was later obtained in [45|46|42|47j . This work also 
pointed out an important connection between the collision kernel in the quantum linear Boltzmann equation 
and a two-point correlation function of the medium known as dynamic structure factor, first introduced by 
van Hove for the description of scattering of a test particle off a macroscopic sample. In the master equation 
the dynamic structure factor, given by the Fourier transform of the time dependent density autocorrelation 
function of the gas, appears operator- valued and accounts for the existence of the canonical stationary 
solution, as discussed in detail in Sect. [8l At the same time, the connection between this quantum master 
equation and the general structure of translation-covariant master equations obtained by Holevo [33j was 
first explored and clarified. 

The final, non-perturbative version of the quantum linear Boltzmann equation, which is discussed in the 
present article, was first derived in [41j . It is noteworthy that this result is not obtained from the weak 
coupling expression by naively replacing the Born approximation of the scattering amplitude with the exact 
expression, since this procedure is not unambiguous. Indeed, the substitution of a^" in (j2.17p by the exact 
scattering amplitude would yield an expression very different from the correct equation as given by (|2.2p 
and (j2.3p . which requires an additional two-dimensional integration. 

Other recent attempts to derive master equations related to the effect of a gas environment are due to 
Alicki [48] and Dodd and Halliwell [49] (see also [50|51l52|53j for recent mathematical work building on the 
covariance under translations). These papers focus on the coUisional decoherence effect for the case of a 
very massive particle, where dissipative effects are not accounted for. In this regard, a reexamination of the 
derivation of the master equation by Gallis and Fleming was presented in [54] . That article discusses several 
ways how to consistently incorporate scattering theory into a Markovian description, which served to ascer- 
tain the numerical value of the localization rate due to coUisional decoherence, thus confirming the prefactor 
implied by [44|46j . The importance of an accurate description of coUisional decoherence was demonstrated 
by experimental measurements with interfering fuUerene molecules j39l55j , which were sensitive to the exact 
saturation value of the localization rate. We note that the techniques used in [54] were also instrumental for 
developing the methods required to combine asymptotic scattering theory with a dynamic description of the 
state evolution, as required in the derivation of the final form of the quantum linear Boltzmann equation. 

Retrospectively, one can distinguish two rather different approaches followed in the various attempts at 
deriving the master equation discussed in the present article. One possible line of argumentation starts 
from the microscopic Hamiltonian for the joint system of test particle and gas, with an interaction term 
specified by a potential depending on the relative distance between gas particle and test particle. A natural 
mathematical framework in this approach is the formalism of second quantization. In this non-relativistic 
quantum field theoretical setting the interaction is naturally expressed as a density-density coupling. In this 
context a second order perturbation expansion leads to the Fourier transform of the density autocorrelation 
function of the gas, that is to its dynamic structure factor [SH]. The result includes statistical corrections, 
which arc naturally accounted for in a field theoretical formalism. However, the restriction to the second order 
term involves the Fourier transform of the interaction potential, corresponding to the Born approximation 
for the scattering cross-section. Going beyond the second order term, thus recovering the full T-matrix 
including statistical corrections, becomes a very difficult task, which has still not been fully accomplished. 
This standpoint has been considered in |57I58I42I49| . 

A second line of reasoning starts from the description of single collision events, formulated by means of 
scattering theory, which has also recently been the object of mathematical investigations |59I60I61I62I63| . 
While this approach admits a non-perturbative description of the effect of a collision, it cannot easily account 
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for degeneracy effects, which may occur if the gas particles are to be considered as being indistinguishable. 
The dynamic effect of the gas is essentially described by weighting the effect of a single collision, obtained 
by tracing over the scattered gas particle, with the rate of collisions. This reasoning was pursued in a 
phenomenological sense in [64l44j . However, since both the rate of collisions and the effect of a single collision 
depend in general on the particle state, it is not straightforward to derive an expression for the infinitesimal 
state evolution which is linear in the density operator, as required in a Markovian master equation, ft is 
achieved in a line of reasoning termed the monitoring approach |65I66| , which is based on concepts drawn 
from the theory of generalized and continuous measurements [28167168] . ft introduces a positive operator 
for the scattering rate and thus yields an expression for the dynamic state evolution, involving the S-matrix 
and square roots of the rate operator, which is manifestly of Lindblad form once the environmental trace 
is done. This approach has been shown to yield nontrivial master equations that have been established by 
other means |65I69| . If applied to the case of a tracer particle in a gas it leads in a stringent calculation to 
the quantum linear Boltzmann equation, as discussed in some detail in |32| . 

Finally, we note that a large number of works can be found in the literature which are related to the 
present subject. While it is impossible to do justice to all authors, we mention a few recent articles re- 
lated to the derivation of a quantum kinetic equation in the spirit of the Boltzmann equation, both linear 
|70l71l7:j|73l74l75l76l77l78l79[Hn] and nonhnear |lllr2ll3ll4ll5ll6ll7l20l81l8-^l83l23l84l85| . 



3. A heuristic derivation 



We now present a heuristic motivation for the particular form of the quantum linear Boltzmann equation 
as given by Eqs. ()2.2p and (|2.3p . It rests upon two basic requirements, its compatibility with (a) the classical 
linear Boltzmann equation, and with (b) the structure of the general form of a covariant completely positive 
master equation. Although our line of reasoning can hardly be called a proper derivation, it is rather 
straightforward and provides an intuitive reasoning why the particular form given above is necessary. 



3.1. The classical linear Boltzmann equation 



It is helpful to first discuss the classical version of the linear Boltzmann equation. In particular, we will 
rewrite it in a form that exhibits a natural correspondence with the quantum version. As already explained 
in Sect. (TJ the classical linear Boltzmann equation describes the collisional dynamics of a test particle 
interacting with a homogeneous, ideal background gas characterized by the momentum distribution ^(p). 
The nontrivial part of the equation is therefore the collisional contribution to the time evolution of the 
distribution function / (P), which we shall denote by 9™"/(P). We start with the form 

i>i°°f(P) = ^j<iP' I V / dP" ("d (P, P) , «1 (P', P')) (3.1) 

where the fundamental role of energy and momentum conservation in each single collision is put in the 
foreground. The differential scattering cross-section ct, which is a function of the relative momenta before 
and after a collision with a gas particle, is taken to satisfy ct (p/,Pi) = cr (Pi,P/), reflecting the invariance 
of the collisional interaction under an inversion of time and parity. We note that the collision integral of the 
original, nonlinear Boltzmann is obtained from ()3.ip by replacing /i (p) with / (p), and by setting M = m. 

As shown in Appendix lA.U by explicitly evaluating the two delta functions in (j3.ip one arrives at the 
expression 

ajo«^,p),r^g:i I ^ I ,ik^^(^rel(k^.P^a)-2 ,el(kj.,Pj.Q) + 3j (3.2) 
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It has already the same form as the right-hand side of Eq. (|2.14p . By exploiting the notation used in Eq. (|2.9p 
and Eq. (|2.12|) one can write (|3.2p more compactly as 

9rV(P) = 1 dQ A/^J (P; Q) / (P - Q) - (P) / (P) , (3.3) 

that is, one arrives at Eq. (|2.13p upon identifying the classical momentum distribution function / (P) 
with (P|p|P), the diagonal matrix elements of the statistical operator in the momentum representation. 
Equation ()3.3p makes it immediately apparent that the classical linear Boltzmann equation takes exactly 
the form of a classical Markovian master equation. 



3.2. General form of a translation- covariant quantum master equation 



The classical result Eq. p.2p already provides a hint on the quantum version of the linear Boltzmann 
equation, in that it tells how the diagonal matrix elements in the momentum representation must look like. 
However, this constraint does obviously not suffice to fix the quantum master equation, which is required to 
describe not only the populations in a given basis, but also the coherences corresponding to the off-diagonal 
elements. 

As is well-known from the theory of open quantum systems [29], the general structure of a Markovian 
master equation is fixed by the Lindblad form, according to which the non-Hamiltonian part must have the 
structure 



(3.4) 



In the present case, however, a more powerful mathematical characterization of the master equation is 
available, which builds on the translational covariance of Cp, and thus provides further cues on the possible 
expressions of the Lindblad operators Lj. Early results on the subject go back to one of the pioneering 
works on quantum dynamical semigroups |86j , and the problem was later reconsidered in |87|88j . The final 
characterization was obtained by Holevo in a series of papers |89|90|91|92j . by building on a quantum 
generalization of the classical Levy-Khintchine formula. 

The linear Boltzmann equation must be invariant under translations, both in the classical and in the quan- 
tum version, because the background gas is taken to be homogeneous, while the interactions are described 
by a two-body potential depending on the relative distance between colliding particles. This property is 
immediately apparent in the structure of the classical equation. As discussed in more detail in Sect. I4.1.H 
it implies on the quantum level that the master equation must be translation-covariant, that is, its action 
must commute with the unitary representation of translations, 

C ( e-'^ P/ Ve^^-'/'') = e-'^■P/'"'£pe'^•P/^ (3.5) 



It was shown by Holevo in the aforementioned work that the structure of a master equation which complies 
with Eq. (|3.5p is given by 



Cp^jdQY. e^'^-^/'^'L, (P; Q) pL] (P; Q) e-'Q'^/" - \ [l] (P; Q) L, (P; Q) , p] 



(3.6) 



with arbitrary functions Lj (P;Q), while an additional Hamiltonian contribution may depend only on P. 
This assumes that the mapping C is bounded and that the underlying statistical process is Poissonian. A 
more general expression including both a Gaussian and a Poisson component has to be considered if one 
also allows for unbounded mappings [47I56|93] . An example of such a Gaussian component is given by the 
master equation Eq. (|5.10p , obtained from the quantum linear Boltzmann equation in the Brownian motion 
limit considered in Sect.O 

Although the form of Eq. (|3.6p complies with the Lindblad result, it provides much richer information 
on the structure of the Lindblad operators. For instance, it states that the only possible dependence on 
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the position operator X can be in the unitary momentum boost operator cxp (iQ • X/h). In particular, the 
diagonal matrix elements of Eq. (j3.6p in the momentum representation are given by 



(3.7) 



(P|£p|P) = / dQ ^ [\L, (Q, P - Q)|' (P - Q|p|P - Q) - \L, (Q, P)!^ (P|p|P) 



3.3. Combining translational invariance with the classical equation 



One can now ask whether a natural choice of the Lindblad operators for the quantum linear Boltzmann 
equation is suggested already by combining the requirements posed by its classical analog and by the charac- 
terization of translation-covariant Markovian dynamics. Specifically, by comparing the diagonal momentum 
matrix elements of the covariant master equation (|3.7[) with the classical expression Eq. p.3p one finds that 
the gain term J^j \^j (Qi P ^ Q)l in Eq. (|3.7[) corresponds to that expression in Eq. (|3.2p which yields the 
classical rate A/f„' (P; Q) in Eq. (j^ . 

It follows that the sum over the discrete set of Lindblad operators (P; Q) must be replaced by an 
integration over a two-dimensional continuous index, in order to recover the five-fold integration in Eq. p.2p . 
Apart from the integral over the momentum transfer Q, this brings about an additional integration over the 
plane of momenta perpendicular to Q, 

J2 Lj (P; Q) pL] (P; Q) ^ y ^ dk^ L (k^, P; Q) pL^ (kx, P; Q) (3.8) 

with L (p, P; Q) the accordingly continuous collection of Lindblad operators. It follows from Eq. (|3.2p that 
they must satisfy 

\L (P, P; Q)!' = (rel (PXQ, Pxq) - §, rel (p^q, Pj.q) + | 

The relation tr (p/,pi) = |/ (p/,Pi)l^ ^'^^d the positivity of p now suggest a natural choice for the Lindblad 
operators, 

/ ^gas p I 



L (p, P; Q) = J^^f rel (p^q, P^q) - ^, rel (p^q, P^q) 




This expression, which involves a complex, operator-valued scattering amplitude, was obtained here in a 
rather straightforward and intuitive way. It is identical with the result (|2.3p . derived stringently by means 
of quantum scattering theory within the monitoring approach. 



3.4. Inequivalent expressions 



The previous section presented a rather suggestive line of reasoning in support of the expression Eq. (|2.3p . 
We emphasize, however, that these arguments are not fully conclusive, since one can formulate quantum 
master equations different from Eq. (|2.2p . which still comply with Eq. (j3.2p and Eq. (|3.6p . Incidentally, 
this is the case for a proposal by Diosi |13], aiming at the generalization of the classical linear Boltzmann 
equation. Clearly, the remaining freedom lies in the fact that the diagonal matrix elements in the momentum 
representation do not uniquely fix the mapping C. The quantum linear Boltzmann equation is an operator 
equation for the density matrix, after all, while its classical version governs only a probability distribution. 
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In order to pinpoint the remaining freedom and the connection to Diosi's result [44j let us return to the 
classical expression Eq. p.2[) . This equation is invariant under a change of integration variables of the form 



k± ^ a(kj., P; Q)ki + 6(kj., P; Q)Pj.q, (3.11) 

with a and b arbitrary scalar functions. This changes the dependence on the test particle momentum P in 
the explicit expression of the classical gain term, such that different master equations are obtained when 
promoting P to an operator and reading off the Lindblad operators, as done above in Eqs. ()3.9p and ()3.10p . 
These various master equations coincide with respect to the diagonal in momentum representation, because 
one can translate back the integration variables when evaluating those matrix elements. 

In particular, by choosing a = I and b ~ m/M the classical linear Boltzmann equation takes the form 



^ ( k_L 



Q JQ 

771 



dk I a ( — ^kj 

m 



Q TO* 

2"' 'ni 



m Q 



M 



-LQ 



TO Q 

toTY 



M 



HQ 



/(P) 



(3.12) 



In this case, the constraint on the Lindblad operators corresponding to Eq. (j3.9[) reads as 



\T I D /-.m2 %asTO /m, Q 771* Q 

P,P;Q =^277^ — P^Q-T' — P^Q + T 
TO^Q \ TO 2 m 2 

/ 7Tl Q TO „\ . 



777, 2 M 

In the same spirit as above, this would lead to the identification 

?^gasTO J, /^TO, _ Q TO, _ ^ Q 



which corresponds to Diosi's result [44| written in operator form. Unlike the quantum linear Boltzmann 
equation (|2.3p . the scattering amplitude is here not operator-valued, since it is independent of P. 

Other choices of a and b can lead to further equations of motion. They all recover the classical linear 
Boltzmann equation, since their diagonal momentum matrix elements coincide, granting in particular the 
existence of a stationary solution. However, they differ in general as far as the description of quantum effects 
is concerned. For instance, the diffusive limit considered in Sect. [5] will lead to different master equations 
of quantum Brownian motion. In this respect, the assignment (|3.10p . and therefore the quantum linear 
Boltzmann equation, is distinguished by the fact that it leads to the minimal correction of the Caldeira- 
Leggett master equation [52], as discussed in Sect. [5] 

As another distinguishing feature of the assignment (|3.10p . the complex rate Mjn (P, P'; Q) resulting from 
Eq. (|2.8p has a physically meaningful form. It is an integration over all those pairs of scattering amplitudes 
/ and /* which are defined by the value of P and of P', respectively, together with the requirement of 
energy and momentum conservation for a given choice of Q and the initial gas particle momentum |66| . 
In an assignment like (|3.14p . in contrast, test particle momenta with different energies 

p2 ^ [P'Y would 

contribute with some common scattering cross-section which does not correspond to the initial and final 
relative momenta determined by P, P', and Q. 
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4. Properties of the quantum linear Boltzmann equation 



We now address features of the quantum linear Boltzmann equation, which are the quantum counterpart 
of properties also exhibited by the classical version. We start by considering the behavior under symmetry 
transformations. Then the relaxation dynamics predicted by the equation is addressed, and in particular 
the approach to an equilibrium state. As in the classical case, the latter can be understood independently 
of the details of scattering cross-section and gas properties. This thermalization behavior is exemplified by 
means of numerical results based on stochastic unravellings which exploit the particular form of the master 
equation. 



4.1. Symmetry transformations 



4.1.1. Spatial translations 

A basic symmetry property of the quantum linear Boltzmann equation is its invariance under spatial trans- 
lations, which was already employed in Sect. [31 One expects this to hold because the coUisional interaction is 
due to a two-body potential which depends on the relative distance between test and gas particle. Together 
with the homogeneity of the gas this implies that a transformation to a frame of reference translated by a 
vector A should not affect the dynamics. This property is expressed at the mathematical level by requiring 
that the mapping yielding the infinitesimal time evolution commutes with the action of the unitary repre- 
sentation of translations, as expressed by Eq. (|3.5p . It implies in particular that if the statistical operator p 
is a solution of the master equation, the translated operator e~^-^'^^^pe^-^'^/'^ also provides a solution. 

To check this invariance of the quantum linear Boltzmann equation let us consider the effect of a translation 
by a vector A on the superoperator Ai appearing in Eq. ()2.ip . We here recall its explicit expression 

TWp = -1 [Ho + i?n (P) ,p]+JdClJ^ dkx [e'Q ></''L (ki, P; Q) pL^ (k^, P; Q) e-^Q'^/'^ 



{p,Lt (ki,P;Q)i(ki,P;Q)} 



(4.1) 



with L (p, P; Q) as in Eq. and Hn (P) as in Eq. ((23)) . 

The mapping in the transformed frame of reference is then given by 



Ma [p] = c'^-^/^M 



-iA-P/h 



(4.2) 



and we have to verify that A4a = -M. 

This is easily checked by separately considering the commutator term and the incoherent part in (|4.ip . 
The Hamiltonian terms are only functions of the momentum operator P and therefore invariant under 
translations. In the incoherent part, the specific form of L (p, P; Q) is not important, since it commutes in 
any case with exp (lA • P/ft), while the unitary operators exp (iQ • X/ft) and exp (— • X/ft) are modified by 
two opposite phases which cancel out. The invariance clearly relies on the special way in which the position 
operator X appears in Eq. (|4.ip . 



4.1.2. Velocity transformations 

The collisional physics underlying the quantum linear Boltzmann equation not only fixes its behavior 
under translations, but also under velocity transformations. This is due to the fact that the interaction rate 
depends on the relative velocity between the tracer particle and the gas particles, the latter being drawn 
from the momentum distribution of the gas particles. A transformation to a reference frame that moves with 
constant velocity V should therefore be equivalent to performing the corresponding shift in the momentum 
distribution of the gas particles. 

Let us therefore consider the superoperator of the infinitesimal time evolution A^V; valid in a frame of 
reference which moves with velocity V compared to the original frame. 
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[p] = c'^^^-^/^M 



-iMV-X/h 



(4.3) 



where M is given by Eq. (|4.ip . According to the microscopic considerations above, we expect that Aiv can 
be obtained as well by evaluating A4 with the transformed gas momentum distribution n (p) jj^ (p) = 
(p — mV). 

To check this symmetry, we first observe that the free kinetic term is left untouched. The explicit form (|2.5p 
of the gas induced energy shift also immediately grants that iJn (P — MY) = (P), where the superscript 
V indicates quantities evaluated with the shifted momentum distribution fjL^ . 

Applying the same analysis to the remaining incoherent part, we note that (|4.3p transforms the second 
term of Eq. (|4.ip into 



>Cv [p] 



dQ / dkj 



^iQ X/h^ (ki, P - A/V; Q) pi' (ki, P - MY- Q) c-'^-^''^ 



\ [p, l' (ki, P - MY- Q) i (ki, P - MY-, Q)} 



The definition of L (p, P; Q) in Eq. imphes 



L (p, P - MY-, Q) = (p + mV^Q, P; Q 



(4.4) 



(4.5) 



Q one 
V In 



By inserting this into Eq. (|4.4p and performing a change of integration variables k^ k'^ = k^ — mV 
also finds for the incoherent part that the transformation (j4.3p is tantamount to the replacement fi fi 
our notation this proves that A^v = , that is, boosting the reference frame is equivalent to translating 
the gas distribution in momentum space. 



4.1.3. Rotations 

On similar grounds one can consider the behavior of the quantum linear Boltzmann equation under spatial 
rotations. Let us denote by R the rotation around a direction n by the angle a. The mapping A4 r in the 
rotated frame is then related to the original mapping M. by 



~mn- J/ h ^^ian - J / fi ^ — ian-J/ h 



(4.6) 



where J is the angular momentum of the test particle, acting as the generator of rotations. 

Usually the collisional interaction between the particles depends only on the relative orientation of in- 
coming and outgoing momenta, / (p/, Pi) = / (Rpf, Rpi)- One expects that the quantum linear Boltzmann 
equation is then invariant under rotations provided the gas momentum distribution is isotropic, i.e., p (p) = 
p (i?p). Indeed, the identity Mb = is then checked easily. 

4.2. Approach to equilibrium 



4.2.1. Detailed balance and the stationary solution 

One of the most important features of both the classical and the quantum linear Boltzmann equation is 
the existence of a stationary solution of the canonical form, e.g. proportional to the thermal momentum 
distribution i^eq(P) given by 



i^EQ (P) 



1 



7r3/2p3 



exp 




(4.7) 



with Pg = y^2M/p. It is reached asymptotically for any initial state and depends only on the mass M 
of the test particle and on the temperature T = 1/ (fee/?) of the gas, provided the latter is described by 
the Maxwell-Boltzmann distribution (|2.6p . We now show that the existence of this solution relies also in 
the quantum case on the fact that the scattering rate M^^ obeys the so-called detailed balance condition. 
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Here and in the following the Hamiltonian part of the master equation does not contribute, so that we can 
concentrate on the mapping £ instead of A4. 

As a first step, consider the quantum linear Boltzmann equation (|2.2[) for a state ;^(P), which is only a 
function of the momentum operator P 

£ HP)] = y dQ [Mf„' (P; Q) u{P - Q) - A/^J (P + Q; Q) u{P)] . (4.8) 

One says that the transition rate A/"^' (P^ P/) = M^^ (P/; P/ — P^) obeys the detailed balance condition 
provided 

p2 \ / p2 \ 



M^' (Pi P2) exp ^-f3^ j = Af^' (P2 -> Pi) exp (^-/^^^ j . (4.9) 

For the gain rate Eq. (|2.9p this implies 

A/^J (P + Q; Q) - A/^; (P; -Q) exp {-f3E (Q, P)) , (4.10) 

with E (Q, P) denoting the energy transferred to the test particle in a collision where the momentum changes 
from P to P + Q, 

i.(Q,P)=(^±^-^ = ^ + ^. (4.11) 

This quantity also relates the Maxwell-Boltzmann distribution at different momenta, 

e-^^(Q^P)i.EQ(P) = i^EQ(P + Q). (4.12) 
Based on the detailed balance condition (|4.9p we therefore have 

>C^^eq(P) = y dQ [Aff„i (P; Q) r.EQ(P - Q) - (P; -Q) z.eq(P + Q)] , (4.13) 

which is equal to zero because the integrand is manifestly odd in Q. If follows that ^'eq(P) is then indeed a 
stationary solution. 

The problem is therefore reduced to checking the validity of the detailed balance condition (|4.9p . This 
can be done considering the explicit expression for the classical rate given by Eq. (|2.10p and taking as the 
stationary distribution function /.j for the gas particles the Maxwell-Boltzmann distribution (|2.6p . Provided 
the scattering cross-section obeys a (p/,Pi) = cr (Pi,P/) one can then directly confirm that ()4.10p holds. 

Finally, we emphasize that the existence of a stationary solution can be proved in a more general and 
direct way by writing the quantum linear Boltzmann equation in terms of the dynamic structure factor of 
the will be discussed in Sect. [51 

4.2.2. Entropy growth 

The celebrated H-theorem for the nonlinear Boltzmann equation [2] states that the function 

H{t) = ldPf{P)logf{P), (4.14) 

satisfies dtH (t) ^ 0. That is, the entropy S [t] = —k-QH{t) associated to the momentum distribution 
function / (P) of the self-interacting gas particles does not decrease under the evolution of the nonlinear 
Boltzmann equation. In particular. dtH {t) = if and only if / (P) ~ f^Q (P), with /eq (P) the canonical 
stationary solution given by the Maxwell-Boltzmann distribution. This tells us that /eq (P) is the only 
stationary solution, reached for i — > 00 independently of the initial condition [94] . 

Unlike the proper Boltzmann equation, the classical linear Boltzmann equation describes an open system, 
implying that the entropy associated to the distribution function / (P) can increase or decrease in general. 
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In this case it is necessary to consider the relative entropy of / (P), defined with respect to another reference 
distribution function g (P), 

i/(/|.9)=y'dP/(P)log^. (4.15) 

As a measure of how effectively one can distinguish the two distribution functions, the relative entropy 
provides information on how close the two distribution functions are. The quantity is non-negative, and 
equal to zero if and only if the two distribution functions coincide. We prove in Sect. lA. 21 that any solution 
/ (P) of the classical linear Boltzmann equation satisfies 

^if(/l/EQ)s$0, (4.16) 

where /eq (P) is the canonical stationary solution given by the Maxwell-Boltzmann distribution. Also in 
this case we have 9tif (/|/eq) = if and only if / (P) = /eq(P). Eq. (|4.16p thus reflects the fact that 
independently of the initial condition any solution approaches the stationary distribution function for t —>■ 
oo. 

The situation is similar in the case of the quantum linear Boltzmann equation. Here one considers the 
quantum correspondence of the relative entropy and studies its behavior under the action of completely 
positive quantum dynamical semigroups [95196197] . 

One should emphasize that the following argumentation can be put on firmer ground only by confining 
the whole system in a region of finite volume. In this case the generator of the dynamics is given by a 
bounded mapping, and the stationary solution of the canonical form is indeed a trace class operator. If the 
test particle can move in an infinite volume, on the other hand, an operator of this form is not of trace class 
since P has a continuous spectrum. The particle is not bounded to a finite region, and it gets completely 
delocalized with elapsing time. 

The quantum analog of the relative entropy is given by the expression 

S {p\w) = fceTrplogp - fceTrplogw, (4.17) 

where p and w are two statistical operators. According to Klein's inequality this quantity is non-negative, 
and equal to zero if and only if the two statistical operators are equal [98|96I97] . 

Taking the quantum relative entropy of the solution of the master equation pt = e*^p with respect to 
the stationary solution peq one is led to consider the time derivative of S (pt Ipeq)- We can now exploit the 
inequality 

S{Up\Uw)i:S{p\w), (4.18) 

valid for any completely positive mapping U [98|96|97] . It provides a quantum analog of the H-theorem for 
dynamical evolutions described by quantum dynamical semigroups. In particular, for h{ = e*^, with Ai 
the generator of a Markovian master equation, and noting that peq is a time translation invariant state, 
e*^PEQ = PeQ; it follows that 

^^S{pt\pKQ)^0. (4.19) 

Eq. (|4.19[) tells us that, even though the variation of the von Neumann entropy associated to the solution 
pt = e*^p does not have a definite sign, reflecting the fact that the system is open, one still has loss of 
information with elapsing time because different initial states get less and less distinguishable. 
In particular, the derivative of the relative entropy with respect to time can be written as 

^5 (ptIpEQ) = fceTr [Mpt (logp* - logpEq)] ■ (4.20) 

It is equal to zero if and only if pt equals the stationary solution, which is unique as a consequence of 
the fact that only multiples of the identity commute with the Lindblad operators exp (iQ -X/fij L (p, P; Q) 
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appearing in the quantum linear Boltzmann equation (j2.ip [95I99| . This result is in complete correspondence 
with the classical case. 

Moreover, according to Eq. (|4.19p any solution of the master equation gets increasingly close to the 
stationary distribution function. This implies in particular that any initial state gets diagonal in momentum 
representation for t oo. That is, the momentum coherences described by the off-diagonal elements decay, 
such that the quantum linear Boltzmann equation goes over to the classical linear Boltzmann equation, 
which describes the remaining populations. 

4.3. Relaxation dynamics 



4.3.1. Momentum and energy relaxation 

The previous section considered the evolution of the statistical operator, which suggested to work in the 
Schrodingcr picture. In contrast to that, we will now use the Heisenbcrg picture when focussing on the 
dynamics of momentum and kinetic energy, that is, on the two observables most relevant for the description 
of a test particle in a gas. 

The Heisenbcrg time evolution of observables requires the adjoint mapping A4*, with respect to the 
quantum linear Boltzmann equation. It is defined through the duality relation 

Tr (AMp) = Tr {pM*A) , (4.21) 

which connects the space of trace class operators, containing the statistical operators p, with its dual, the 
space of bounded observables A. The dynamics of a Heisenbcrg observable At is thus given by 

^At = 4:[At,Ho + i?„(P)]+rA, (4.22) 
dr iri 



with 



£*A = y" dQ y dk^ (k^, P; Q) c-'^-^/r^p^^^Q X/hj^ ^^^^ p. 
-^{l' (ki,P;Q)L(kx,P;Q),A} 



(4.23) 



We will take the Heisenbcrg operator to coincide with the corresponding Schrodinger observable at t = 0, 
i.e., Ao = A. 

An important class of observables are those which are only functions of the momentum operator. At = 
At (P). The equation of motion in the Heisenbcrg picture then simplifies considerably, 

^At = y dQ dki \L (kx, P; Q)|' [c-'Q-^/''Atc*Q-^/'"' - At] . (4.24) 

Note that C, and equivalently £*, is covariant under translations in the sense of Eq. p.Sp . which implies 
that the algebra generated by the momentum operator is left invariant. Consequently, Aq ~ Aq (P) implies 
At = ^t(P) for i > and therefore [At, P] = 0. Recalling the definition of the quantum and the classical 
gain rates in Eqs. (|2.8p and (|2.9p . we find that observables given by a function of momentum obey 

^^t(P)= y dQK;(P + Q;Q)[^(P + Q)-At(P)], (4.25) 

in strict analogy with the classical formulation. 

Let us now focus on the time evolution of the expectation values of momentum and kinetic energy, 
A (P) = P and A (P) = P^/ (2Af ). Using the momentum representation one can evaluate Eq. (|4.25[) for these 
quantities by working with real numbers rather than operators. In particular, for A(P) = P the quantity 
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[A(P + Q) — j4(P)] appearing at the right-hand side of Eq. (|4.25p is simply the momentum transfer Q in 
the sing le colhsion. For A{P) = P^/ (2A/) it is the energy transfer i?(Q,P) as given by Eq. (|4.1ip . The 
positivity of M^^ thus allows one straightforwardly to confirm from Eq. (j4.25p that the change of momentum 
is positive when the momentum transfer is positive. Likewise, a positive energy transfer increases the energy. 

The dynamic equation (|4.25p for the momentum and kinetic energy can be evaluated explicitly if one 
assumes a constant scattering cross-section |/(p/,Pi) P = o'tot/47r and a free gas of Maxwell-Boltzmann 
particles. Inserting the expression of Afj^' into Eq. (|4.25p it can be shown in this case [lOOj that the time 
evolution of the momentum expectation value is given by 

^(P). - -A'-^r.iP.F, ,4.26) 



Accordingly, the kinetic energy E = P^/ (2Af) evolves as 

d 16 m, , ^ / 1 5 „^ m \ ^ 3 ^ / 3 3 ^^m\. „„n 

d? - - M ^^^^^^ (- 2 ' r ~^^M ) ^ - 2j3^^^^ [-r r -^^M ) ^'-^'^ 

In these equations we have denoted the most probable velocity of the gas particles as vp ~ pp/m = 
yj2/ (/3m). Moreover, 

T/j = ngasW/30-tot (4.28) 

is the total scattering rate for a flux of particles with the most probable velocity w^, and the symbol li^i 
denotes confluent hypergeometric functions. Their explicit expressions are given by 

and 



I -.x^) = 1 I [5 + 2x-] e-= + [3 + 12.x^ + 4x^] , (4.30) 

which are both positive monotone increasing functions. From these expressions it is immediately clear 
that in general there is no closed evolution equation for either the first or the second moment of the 
momentum operator, (P)pt or (P^)^^, since moments of arbitrary high order are involved in the equation, 
due to the presence of the confluent hypergeometric functions. This becomes manifest in strong deviations 
of the moments from Gaussian statistics, as can be observed by studying the time dependence of high order 
cumulants of the momentum distribution |101| . 

However, the equations do get closed in the limit of a very massive test particle close to thermal equilibrium, 
which corresponds to the diffusive limit considered in Sect.[5l In this case, the velocity V ^ P/M of the test 
particle is taken to be much smaller than the typical velocity vp of the gas particles. Thus taking V/vp <^ 1 
and m/Af <C 1 in Eq. (|4.26p and Eq. (|4.27p the confluent hypergeometric functions can be replaced by unity, 
since iFi [a, 7; 0) = 1, while the reduced mass to* is replaced by m. This leads to 

^(P)p.=-^(P)p* (4.31) 

which describes a friction proportional to velocity, leading to the expected exponential relaxation to a 
vanishing mean momentum. Similarly, one finds 



A(E)„.-2„((E)„-1 



E>„=-2i( (E)„- — ) (4.32) 



which shows that the mean kinetic energy relaxes exponentially to the equipartition value 3/2kBT. The 
friction coefficient ?/ in Eq. (|4.3ip and (|4.32p is given by 
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(4.33) 



4.3.2. Simulation algorithm 

The complexity of the quantum Uncar Boltzmann equation prevents its solutions from being analytically 
tractable, even if one allows for special choices of the scattering amplitude and the initial state. Nevertheless, 
one can efficiently simulate the temporal behavior of various interesting quantities by using quantum trajec- 
tory methods, which build on the Lindblad structure of the master equation and its translation-covariance. 
These techniques are based on stochastic differential equations for the state vector \tp (t)), whose solutions are 
a stochastic process in the Hilbert space. A stochastic unravelling of the quantum master equation is achieved 
if the ensemble average E {t)){ip (t) |] yields a solution of the master equation [1021103110411051^ . 

One way of motivating these unravelling techniques is to consider measurement schemes, which lead to 
instantaneous changes of the wave vector conditioned on the random outcomes. Together with the covari- 
ant structure of the master equation described above, this consideration leads to a natural choice of the 
unravelling. It is a piccewise deterministic process, described by a deterministic time evolution, which gets 
interrupted by random jumps which change the particle momentum. This unravelling corresponds to a 
measurement scheme in which one monitors the momentum transferred between system and environment. 

For the sake of simplicity we will restrict our analysis to the Born approximation of the equation, where 
the scattering cross-section depends only on the momentum transfer. Although this is not a conceptual 
restriction, it leads to major simplifications in the sampling of the jump events. 

In Born approximation the quantum linear Boltzmann equation reads 



±p = -'-[H,p]+ J dQ e'Q>^/''LB (P; Q) pL^ (P; Q) e-'Q'^/" - ^ (P; Q) Lb (P; Q) , p} 



(4.34) 



where we have used the more compact notation H = Ho + i^n (P) for the Hamiltonian term, and the Lindblad 
operators are given by Eq. (|2.17p . 

In the Monte Carlo wave function method [105|29j one then considers a stochastic differential equation 
for the stochastic wave vector tp (i), which has the form 



dli^it)). 



h^^^^\j dQI|iB(P;Q)|V'(t)) 

^.Q.X/r.^^ (P; Q) 1^ {t)) 



dQ 



||Ls(P;Q)|V'(0>ll 
Here, the rate operator G is given by 

G = j dQ4(P;Q)LB(P;Q). 

The field of Poisson increments dNq (t) satisfies 

dNq {t) dNq, (t) = 5^ (Q - Q') dNq (<) 
E[dNq {t)]^\\LB{P;Q)\^m?^t. 



IV' (0) 



\^:{t))dt 
dNq {t) . 



(4.35) 



(4.36) 



(4.37) 



where the symbol E denotes the ensemble mean over realizations of the process. The solutions of (|4.35|) 
starting from an initial state \'ip (0)) provide an unravelling of the master equation (j4.34p in the sense that 



Pt = E[|V;(t))(^ {t) 



(4.38) 



ifpo = |^(0))(V'(0)|. 

In spite of the complicated form of the stochastic differential equation, which involves a field of Poisson 
increments, the Monte Carlo wave function method provides a rather simple simulation algorithm. The 
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realizations of the process are characterized by deterministic time evolutions interrupted by jumps. If a 
jump into the state {ip (t)) occurs at time t, it will evolve until the time t + t of the next jump according to 
the deterministic equation 

The waiting time r between two jumps is determined by the cumulative distribution function 

F{t) = 1- ||e-*^^^/'*-^^/2|^/. {t))f. (4.40) 

A straightforward way to obtain the random variable t numerically is to draw a random number rj from a 
uniform distribution over the interval [0, 1] and to invert the relation \\c^^^'^^^^^'^^^\tp {t))\\'^ = rj. The jump 
is then implemented by the replacement 

where the unitary operator exp (iQ • X/fi,) effects a momentum kick of amount Q. These momentum transfers 
must be drawn from the probability distribution 

which is normalized by the total scattering rate out of the state ^ (t) , 

(t) \G\yj it)) = J dQllLs (P; Q) IV it)) (4.43) 

The random jumps in this unravelling are therefore naturally associated to the collisions with the gas 
molecules experienced by the test particle. The realizations of the process are thus given by a deterministic 
evolution interrupted by random momentum kicks corresponding to collisions. The value of the momentum 
transfers and the rate of collisions arc fixed by the state dependent quantities Eq. (|4.42p and Eq. (|4.43p 
respectively. 

The described algorithm can be implemented in a particularly simple way if the initial state is a discrete 
superposition of momentum eigenstates, |'0(O)) = c^ lPj) . Due to the translation-covariance (|4.2p of 

the master equation superpositions of this kind are preserved during the stochastic time evolution in the 
sense that the quantum trajectory remains a superposition of N momentum eigenstates at all times. The 
stochastic process is thus mapped into the space of the different amplitudes Cj and momenta Pj . 



4.3.3. The loss term 



When discussing the momentum and energy relaxation in Sect. 14.31 it was not necessary to work out the 
loss term of the quantum linear Boltzmann equation explicitly, since it cancels out in Eq. (|4.25p . However, 
this expression is of independent interest; it appears in the Monte Carlo simulation of the quantum linear 
Boltzmann equation through the rate Eq. (|4.36p , and it also determines the loss of visibility due to coUisional 
decoherence in matter wave interference experiments, as we shall discuss in Sect. [71 We therefore provide an 
explicit expression assuming that the cross-section depends only on the relative velocity, which includes the 
case of a constant cross-section discussed in Sect. 14.31 

The classical loss term in the quantum linear Boltzmann equation is given by 

AossP=-i{A/ot(P),p}, (4.44) 

as evident from Eq. (|2.7p . It is convenient to rewrite the classical loss rate (P), in terms of an effective 
scattering cross-section aes (P) defined as the area to be multiplied by the current density rigi^gV of the test 
particle. 
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p 



(4.45) 



The effective cross-section can be calculated exactly if the scattering cross-section is independent of the 
scattering angle (s-wave scattering) , and if it depends on the energy only as a power of the relative velocity 
between test particle and gas particles. The microscopic cross-section can then be expressed as 



o- ( rel (p±q,P±q) - y,rcl(piQ,P_LQ) + y 



This yields the effective scattering cross-section |106| 



(Teff (P) = SV^T 
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(4.46) 



(4.47) 



where i-Fi denotes the confluent hypergeometric function and = ■\/2/ (/3?7i) the most probable velocity 
in the gas. 

For a = one recovers the case of a constant scattering cross-section \f (p/, Pi)|^ = crtot/47r, used below to 
perform a Monte Carlo simulation of the quantum linear Boltzmann equation. The corresponding expression 
is 



(Toff (P) = crtot—;=—p:-iJ'i 



hi P 




where 
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[1 -f 



/TT erf (x) 



For small velocities of the test particle V jv^ ^ \ the expression Eq. (|4.48p simplifies to 

v^vp 2 Mvfi V 

Ccff > atot—F= T^^^ftotTT, 

Vtt P V 



(4.48) 



(4.49) 



(4.50) 



with V = -y/S/Trm/? the average thermal velocity. This corresponds to the famous 1/V dependence of the 
total scattering cross-section for small velocities, expected when the scattering probability is independent of 
the velocity of the tracer particle, as verified in neutron scattering [7]. 



4.3.4. Simulating the momentum and energy dynamics 

We now discuss a numerical study of the relaxation properties of the quantum linear Boltzmann equation 
[101| . which is based on the unravelling given by Eq. (|4.35p and the algorithm presented above. The use of 
dimensionless variables, naturally suggested by the typical physical lengths of the problem, turns out to be 
very convenient. Specifically, we introduce the variables 



and 



U 



K 



V 



Q 



(4.51) 



(4.52) 



They provide the momentum of the test particle, scaled with respect to the most probable velocity Vfj in 
the well as the scaled momentum transfer in a single collision. 

For the sake of simplicity we also focus on the case of a constant cross-section. In terms of the new variables 
the Lindblad operators then read as 
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Fig. 4.1. (Color online) Single quantum trajeetory realizations of the scaled momentum U(t), as described in Sect. 1473.41 The 
panels differ in the ratio m/M between the gas and particle mass, which decreases from top left to right bottom (top row: 
m/M = 2 and m/M = 1, bottom row: m/M = 0.5 and m/M = 0.1). One observes that the number of collisions required to 
get into the thermal cloud increases as m/M gets smaller, while the width of the cloud decreases. 



iB(U;K) 
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/3 



exp 



K 

T 



U • K' 

K , 



(4.53) 



We assume that the initial state is a momentum eigenstate |U). labeled by the dimensionless momentum 
(|4.5ip of the test particle. The algorithm considered in Sect. 14.3.21 then leads to a pure jump process in U, 
whose realizations can be obtained by means of the standard algorithm used for the stochastic simulation of 
classical Markovian master equations [107] . This is due to the fact that the deterministic evolution given by 
Eq. (|4.39p only affects the normalization in this case, since both H and G are only functions of the momentum 
operator. 

According to Eq. (|4.40p and Eq. (|4.43p the waiting time distribution for a jump is characterized by the 



cumulative distribution function F (t) = 1 — exp — F {U) r 



in our case, where T ([/) 



Mf^, (U) is the 



dimensionless version of the loss rate. According to Eq. (|4.45p and Eq. (|4.48p it reads 
2 / 1 3 



T{U) 



- 13- 



--iFi 



2 2' 



(4.54) 



with as in Eq. ((428)) and iFi as in Eq. (|449)) . 

Each jump is characterized by a momentum transfer drawn from the probability density (|4.42p . In the 
present case, this distribution depends only on the modulus K of the momentum transfer and on the cosine 
^ = U • K/ {UK) of the angle between U and K, 



2V^r(c/) 



K eyij) 



(4.55) 



while components of the momentum transfer perpendicular to U are uniformly distributed. For each jump 
one has to draw random numbers (A', ^) according to the probability density (|4.55p and perform the shift 
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Fig. 4.2. (Color online) Averages of squared mean momentum (top) and mean squared momentum (bottom) for the mass 
ratios m/M = 1 (left) and m/M = 0.2 (right). The solid line gives the results obtained from the Monte Carlo method, 
while the dashed lines correspond to the analytic predictions in Eq. 14.5711 and Eq. 14.581 1. respectively. As one expects, the 
analytic approximation deteriorates if one increases the mass ratio (left panels). In this regime the decay rates for the (squared) 
momentum and the kinetic energy also no longer coincide. 



U ^ U + Km* /A/ with K given by 



U 



U X e 



^ ^ U X e 



(4.56) 



with e a random unit vector. 

Single reaUzations of the process for different values of the mass ratio are shown in Fig. 14.11 as three- 
dimensional plots. One observes that the ratio between the masses of the gas and test particle affects the 
number of collisions which are necessary in order to drive the test particle momentum close to its equilibrium 
value. At smaller mass ratios m/M, an increasingly large number of collisions is required to thermalize the 
test particle. The fluctuations around the equilibrium value, on the other hand, decrease in this limit. 

By averaging over a sufRciently large number of realizations one can study the relaxation of the momentum 
and of the kinetic energy of the test particle. According to the scaling introduced in Eq. (|4.5ip the expected 
equilibrium value 3/ (2/3) of the kinetic energy corresponds to (U^)eq = 3m/ (2M), while the momentum 
should relax to (U)cq — 0. This behavior is confirmed well by the numerical simulations reported in Fig. 
14.21 where the squared mean momentum (U)gq and the mean squared momentum (U^)oq arc plotted as a 
function of time in units of the inverse of the total scattering rate F/j. 

The numerical results in Fig. 14.21 are compared to the analytic predictions Eq. (|4.3ip and Eq. (|4.32p for 
the mean momentum and mean kinetic energy respectively, which are valid for m/M <C 1. In terms of the 
dimensionless variable U the solutions of Eq. (|4.31[) and Eq. (|4.32[) read as 



(U)?,=(U)^„e 



(4.57) 



and 
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(4.58) 



with T] given by Eq. ()4.33p . As one observes from Fig. l4.2l the analytic approximation is indeed very good for 
small mass ratios m/M <^ 1, while deviations appear for higher values of the ratio. These changes show up 
in particular through the appearance of different decay rates for the momentum and the kinetic energy. For 
large m/M one can also observe strong deviations from Gaussian statistics, e.g., by studying the behavior 
of cumulants of order higher than two [lOlj . 

5. The limit of quantum Brownian motion 

One of the most important limits of the quantum linear Boltzmann equation is the situation where the 
test particle has a much greater mass than the gas particles. In the classical case, this limit turns the linear 
Boltzmann equation into the Kramers equation for the description of Brownian motion. As discussed in the 
following, the situation is quite similar in the quantum case, where one obtains a quantum Fokker-Planck 
equation describing quantum Brownian motion. 

The term quantum Brownian motion is generally used for master equations intended to phcnomcnolog- 
ically describe the quantum analog of classical Brownian motion [108110911101111] . without accounting for 
the microscopic details of the underlying dynamics. In contrast to such models, which usually assume a 
linear coupling to a harmonic bath, we focus in the following on the specific situation of a massive test 
particle suffering collisions with a surrounding Maxwell-Boltzmann gas. The microscopic expressions of the 
friction and diffusion coefficients arc worked out as a function of the cross-section, and we put into evidence 
the new quantum terms appearing in the master equation corresponding to the Kramers equation. 

5.1. The classical Kramers equation 

In order to discuss the Brownian motion limit of the quantum linear Boltzmann equation it is useful to 
first recall the relationship between the classical linear Boltzmann equation and Kramers' Fokker-Planck 
equation. It was Rayleigh who first investigated the limiting case of a very heavy test particle in the collision 
integral of the linear Boltzmann equation |112|113j . The resulting equation in this so-called Rayleigh limit, 
or Brownian limit, is the well-known Kramers equation. For the homogeneous case and in the absence of an 
external force it reads as 



with (7(i^; upp) the differential scattering cross-section as a function of the scattering angle and of the 
momentum of the gas particle p = upp. 

The derivation of the microscopic expression of the friction coefficient is not easy to find in the literature. 
We therefore briefly recall the derivation of (|5.ip and (|5.2p in Sect. IA.3I for the reader's convenience and 
for the sake of comparison with the quantum calculation. The crucial assumption required to derive the 
Fokker-Planck Eq. (|5.ip from the classical linear Boltzmann equation (|3.ip is that the test particle is much 
heavier than the gas particles, m/M ^ 1, and that its momentum is close to the typical thermal value 
Pf3 = \/2M/p. The collisions with the gas particles are then characterized by small momentum and energy 
transfers. The typical value of the momentum transfer Q is of the order of the momentum of the colliding 
gas particles = y/2m/ f3, and therefore much smaller than the momentum Pp of the test particle, while 
the energy transfer is on the order of E (Q,P) ~ 1//3 i/m/M, which is smaller than the thermal energy 1//?. 
As discussed in Sect. IA.3[ under these conditions one can expand / (P') in Eq. (|3.ip around P for small 
momentum transfers Q = P' — P, and arrive at Eq. (j5.ip by sorting the orders in an m/M expansion. 




(5.1) 



The friction coefficient 77 is given by [1141115] 




di9 sini? (1 — cosi}) a{-d; upp), 



(5.2) 
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5.2. The quantum Brownian limit 



5.2.1. Conditions of validity 

The situation is more complicated in the quantum case, since one is dealing with operators whose values 
can be estimated in a meaningful way only by considering suitable matrix elements. We will show that 
the quantum master equation which corresponds to the classical Fokker-Planck equation is obtained by 
expanding the scattering kernel up to second order in the operators X and P. It is important to stress that 
the result is not equivalent to a direct application of the correspondence principle on the classical result, 
which would lead to a master equation which does not generally preserve the positivity of the statistical 
operator |116I117| . 

The expansion can be performed provided the test particle is very massive and not far from thermal 
equilibrium, that is to say, close to diagonal in momentum representation. More precisely, the off-diagonal 
elements (P|p|P') may differ significantly from zero only for AP := |P - P'l < ^/M/p. In the position 
representation, the validity of the expansion requires that the motional state exhibits coherences only over 
a length scale of the order of the thermal wavelength of the test particle, such that (X|p|X') is appreciably 
different from zero only within a range given by the thermal de Broglie wavelength Ath = \/2TTh?l3/M ^ i.e., 
for AX := |X — X'| < Ath- The small mass ratio m/M <^ 1 implies in particular that when the test particle 
gets close to thermal equilibrium its coherence length of the order of its thermal de Broglie wavelength is 
much smaller than the coherence length of the environment given by the thermal wavelength of the gas. 



5.2.2. Diffusive limit of the quantum linear Boltzmann equation 

In order to study the Brownian motion limit of the quantum linear Boltzmann equation, we start with the 
incoherent part given by Eq. (|2.2p . Using the explicit expression of the Lindblad operators Eq. (|2.3p with 
the distribution function (|2.6p of a Maxwell-Boltzmann gas, one obtains the exact expression 
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(5.3) 



In the expansion up to second order in P the contributions from the terms exp /JmP^g/ (4A/^) j cancel out, 
such that they can be replaced by unity. We further expand the various quantities under the assumption that 
the change in momentum of the Brownian particle is small compared to the scales involved in its motional 
state. More specifically, one has to assume that for typical values of the momentum transfer Q the relevant 
matrix elements of the statistical operator vanish unless 



-AX<1 



(5.4) 



and 



PQ 

M 



AP<1. 



(5.5) 



These conditions are both satisfied if the Brownian state is close to thermal and M 3> m, such that to* « to, 
since the momentum transfer Q is then typically of the order of the momentum = ^/2/rn/P of the colliding 
gas particles. This implies in particular that 
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and 



M 



(5.7) 



Apart from the exponentials in Eq. (|5.3p one also has to expand the scattering amplitude, which appears 
as an operator-valued expression. However, for the assumed small mass ratio the argument of the scattering 
amplitude becomes 



rel(k^,P^Q)-|«k^-^P^Q-§, 



(5.8) 



such that the leading term in an expansion in m/M does not depend any longer on P_lq. This removes an 
operator dependence arising from the scattering amplitude. By further expanding the exponentials involving 
Q • P in Eq. ()5.3|1 one is then left with 
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(5.9) 



We here omitted the term linear in the momentum transfer, since it vanishes upon integration due to 
the isotropy of the integrand. By exploiting also the invariance of the scattering cross-section under parity 
transformations the quantity QjQk can be replaced by Sji^Q^ /3. We have thus put into evidence the operator 
structure of the limiting expression of the quantum linear Boltzmann equation for the case of quantum 
Brownian motion. 
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(5.10) 



This equation is still covariant under translations |47I56I93] . since the limiting procedure does not spoil the 
symmetry properties. 



5.2.3. Microscopic expression for the friction constant 

We now evaluate the friction coefficient 77 in Eq. (|5.10|) . which is defined by 
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(5.11) 



The exact expression of 77 will turn out to be the same as the classical friction coefRcient (j5.2p . Comparing the 
procedure followed here with the calculations in Appendix IA.3I it turns out that the microscopic derivation 
of the friction coefRcient is more straightforward in the quantum case. Upon using (5 (p • Q) =6 (p||) /Q and 
introducing the new variables = p+Q/2 and p/ = p — Q/2 the expression for rj becomes 



77: 
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(5.12) 



As a final step, we introduce the angle d between pi and p/ as an integration variable, and recall that 
the scattering cross-section is assumed to depend only on pi and thus finally obtaining for a Maxwell- 
Boltzmann gas distribution 
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with u as in (|5.2p . This expression coincides with the classical result |114| . recalled in Eq. (|5.2p . For the case 
of a constant scattering cross-section it reads as 

which also agrees with the result (|4.33p obtained in Sect. 14.31 

5.3. The master equation of quantum Brownian motion 

Putting the master equation for quantum Brownian motion Eq. (|5.10p in the standard form, it finally 
reads 

i=i i=i j=i 



The friction coefficient 77 is given by Eq. (j5.13p . while Dpp has the meaning of a diffusion coefficient 

M 
1 

and Dxx is related to Dpp and to 77 by 



Dpp^V—, (5.16) 



This equation admits an exact solution, though the explicit expression is a bit cumbersome }118| . We shall 
consider it in Sect. 17. 4( when studying how this master equation describes decohcrence effects. 

We note that the Caldeira-Leggett master equation, describing the high-temperature regime of the Caldeira- 
Leggett model |119I120I121I111| . is related to Eq. (|5.15p by neglecting the third term in (|5.15p . i.e., by setting 
Dxx = 0. This equation, which cannot be brought into Lindblad form, is also obtained from the classical 
Kramers equation if one takes the naive operator correspondence d/dPj Xj / (ih). 



5.3.1. The diffusion coefficients and complete positivity 

The diffusion coefficient Dpp can also be related to the second moment of the thermal momentum {P'^)i3 
3M/ j3 according to 



Dpp = \{P^)p, (5.18) 

while Dxx, which is sometimes denoted as "quantum diffusion coefficient" or "momentum localization coeffi- 
cient" , can be written by means of the thermal de Broglie wavelength Ath = \J2ttW(3/M of the test particle 
as 

^x. = r^A?h. (5.19) 

As discussed in |122I123I124I125I126I29| a necessary and sufficient condition for master equations of the 
form (|5.15p to yield a completely positive dynamics is 



DxxDpp-^-!^, (5.20) 



16 

since only then they can be cast into Lindblad form. This is an important statement, because in many 
approaches the coefficients t], Dxx, and Dpp are treated as free, phenomenological parameters. In contrast 
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to that, the derivation from the quantum hncar Boltzmann equation yielded microscopically formulated 
definitions and relations of the coefficients (|5.15p - (|5.17p . It is therefore reassuring that D^x and Dpp from 
Eq. (|5.17p and Eq. (|5.16p satisfy Eq. (|5.20p with an equality sign. 

For a given choice of 77 and Dpp = rjM/P the coefficient D^x obtained in Eq. (|5.17p therefore takes the 
minimal value allowed within the framework of a Lindblad master equation. This implies that the equation of 
quantum Brownian motion obtained from the quantum linear Boltzmann equation is as close to the classical 
Kramers equation as it can possibly be, within the bounds set by the requirement of complete positivity. 

5.3.2. The quantum diffusion term 

It was discussed by several authors |122I123I124I125I127I29| that the lack of complete positivity of the 
Caldeira Leggett master equation can be corrected by adding a "position diffusion term" of the form 
[Pj, [Pj,p]]. In view of the fact that the Brownian limit of the quantum linear Boltzmann equation yields 
the minimal correction, it is worth noting that this is not achieved in other microscopic approaches. 

Diosi's proposal [U], for instance, leads to an expression of the coefficient Dxx, whose relation to Dpp 
is not given exclusively by thermodynamic quantities. Rather, it involves a complicated average over the 
scattering cross-section. More refined derivations within the Caldeira Leggett model, on the other hand, 
which seek for correction terms beyond the high temperature approximation |128|129j . lead to a coefficient 
Dxx which exceeds the present, minimal value by a factor of 4/3. 

As already mentioned, the term J2j=i i^j^ l^j^P]] ^'i- l|5.15p has no counterpart in the classical Fokker- 
Planck equation (|5.ip . It describes a position diffusion process, which is expected to appear together with 
momentum diffusion, as implied by the inequality Eq. (j5.20p . That this is a pure quantum phenomenon is 
confirmed by the fact that the coefficient Dxx depends explicitly on h; it vanishes in the semiclassical limit 
0. 

Even though the occurrence of this quantum effect leading to momentum localization was predicted in 
many theoretical papers, it is still unclear how to observe the effect experimentally. The main difficulty lies 
in the fact that the position localization term, with its direct classical analog, dominates under typical ex- 
perimental conditions, due to its strong temperature dependence, while the quantum diffusion term vanishes 
in the limit of high temperatures. Effects of the latter might show up by either inducing decohcrence with 
respect to the momentum basis, or by looking for corrections to the classical diffusion coefficient. The last 
possibility was considered for the strong friction limit of the quantum master equation |130j . Compared to 
the classical case, the Einstein diffusion coefficient appears multiplied by a factor 1 -I- (r](3h) /16, involving 
the ratio between two characteristic times (3h for the bath of gas particles and 1 /rj for the system. The sep- 
aration of time scales typical for a Markov description of the dynamics implies that this correction should 
be extremely small under typical conditions. 

6. Matter wave optics in gaseous samples 

So far, we mostly dealt with features of the quantum linear Boltzmann equation that are also shared by 
its classical counterpart. In particular, its symmetry and relaxation properties, as well as the existence of 
the stationary solution, are already described by the behavior of the momentum populations, i.e., of the 
diagonal matrix elements of Eq. (|2.7p . In the following two sections, we will now consider properties of the 
quantum linear Boltzmann equation that become manifest in interference phenomena, where it is decisive 
that the motion of the test particle may be in a coherently delocalized, wavelikc state. 

An important quantum effect that one expects to be described by Eq. (j2.2p is the dynamical transition 
of an initially delocalized superposition state into an incoherent, particlelike mixture. Such decohcrence 
phenomena are described in Sect. [T] But before that, we focus on the dynamics described by the coherent 
part of the quantum linear Boltzmann equation and its loss term. They play an important role in experiments 
involving spatially separate beams. 

The claim that massive particles can exhibit wave- like behavior is one of the central predictions of quantum 
mechanics; it is confirmed experimentally with neutrons [1311132] . atoms |133ll34l35j . and even complex 
molecules [135j . In such experiments, the de Broglie wave length AdB = h/ (MV) of the matter waves is 
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usually much shorter than the scale of change of the potential produced by the gratings and coUimation slits, 
which act as optical elements. It is then often permissible to neglect diffraction phenomena by describing the 
matter wave in terms of well defined beams, which are coherently split and rejoined by means of interference 
gratings. The effect of an additional external potential on those beams can then be accounted for, in complete 
analogy to the optical case, by introducing an index of refraction. 

In a number of Mach-Zehnder-type interference experiments a truly macroscopic separation of the inter- 
fering beams was achieved. This makes it possible to perform experiments where the interfering particles 
interact with an external system in just one arm of the interferometer. The modification of the interference 
fringe pattern then provides sensitive information of the strength of interaction. Its theoretical description 
requires the use of quantum kinetic equations if the interaction is due to a dilute background gas, as realized 
in experiments with sodium and lithium atoms [37|136|38] , or due to a homogeneous sample of thermalized 
condensed matter, as in the experiments performed with the perfect crystal neutron interferometer [1371132] . 



6.1. The gas induced energy shift 



The presence of a background gas manifests itself not only in the incoherent part C of the quantum 
linear Boltzmann equation, but also in a modified Hamiltonian H = P^/ (2M) + Hn (P) appearing in its 
coherent part. The additional term Hn (P), which has been neglected so far, describes the modification of 
the free matter wave dispersion relation due to the coherent interaction with the background gas. It becomes 
observable in the abovementioned Mach-Zehnder interferometers as a gas- induced phase shift. 

According to Eq. (|2.5|) . derived in [32 j the energy shift (P) is related to a "forward scattering process" 
which does not change the momenta of the gas and the test particle. The forward scattering amplitude 
/o (p) = / (p, p) appears in a thermal average over the momentum distribution of the gas, 

(/o (P))^ = 1 dp flip) f (rel (p, P) , rel (p, P)) , (6.1) 

allowing us to write more compactly 

Hn{P) = -27rh^'^Re{foiP))^. (6.2) 

To observe the effect of Eq. (j6.2p an experimental situation is required where the interfering matter wave 
beams are strongly collimated and directed into narrow apertures. In this case any proper, momentum 
changing collision due to a dilute background gas will unavoidably lead to a loss of the particle. This implies 
that all particles arriving in the detector experienced at most the effect of "forward scattering" . It can be 
observed as a phase shift in the interference pattern, as done in the experiments [371136138] . 

If the particles do not get lost, their motional state is modified at most coherently in the described experi- 
mental situation. The purity of the matter waves impinging at the detector therefore remains unchanged by 
the presence of the background gas, rendering a description by means of a wave equation possible. The effect 
of particle loss can be accounted for by an attenuation of the amplitude, as argued below in Sect. Inci- 
dentally, the loss term introduced in Eq. (|2.13p and discussed in Sec. 14.3.^ can be related to the imaginary 
part of Eq. (|6.ip . This fact naturally leads to the introduction of an optical potential. 



6.2. The optical potential 



It is helpful to express the damping of the wave amplitude due to collisional losses in term of the forward 
scattering amplitude. Let us start from Eq. (|3.1|1 for the loss term in the classical linear Boltzmann equation. 
Switching to center-of-mass variables with respect to p' and P', denoting k = rel(p',P') and K = p' + P' 
one obtains 
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(P) = !^ J dkk J dp/^(p)(5(fc- |rel(p,P)|) J dfikcr(k,rel(p,P)). (6.3) 
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We also introduce the total scattering cross-section 

CTtot (p) = J df^n cr {pn, p) , (6.4) 

where n is a unit vector with df2n the associated solid angle element. By exploiting the optical theorem 
[138], 

4:Trh 

ctot (p) [/ (P' P)] ' (6-5) 

and by using the thermal average defined in Eq. (j6.ip , one finds that the loss term takes the simple form 

M^l, (P) = ^47rHm(/o (P)>^. (6.6) 
Recalling Eq. (|2.7[) and Eq. (|6.2p the quantum linear Boltzmann equation can therefore be written as 



A - i 
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^Anmm{,k{P))^,p} (6.7) 



JdQj dk^ e^Q-^/'^'L (ki, P; Q) pit (k^, P; Q) e-''^' - 



This suggests to introduce a non-hermitian operator, often denoted as optical potential, of the form 



Vopt{P) = -27rh'^{fo{P))^. (6.8) 



One thus arrives at a compact and suggestive expression for the quantum linear Boltzmann equation 
d_ 

dt^ ih 



% = 1 [Ho, p] + \ (Kpt (p) p - pkU (p; 



+ J'^qJ dki c^'^-^/'^L (ki, P; Q) pLt (k^, P; Q) e-'^■^/^ (6.9) 

The first three terms at the right-hand side describe a coherent dynamics. They preserve the pure states, 
but not their normalization, due to the fact that the optical potential is not self-adjoint. The completely 
positive mapping expressed by the last term describes an incoherent contribution. It transforms pure states 
into mixtures and accounts for the preservation of the trace. 

6.3. Index of refraction 

As described above, typical matter wave interference experiments are characterized by a well-collimated, 
stationary beam of particles, with a momentum centered around a reference value P = HK. in the longitudinal 
direction. One can then exploit the formal analogy between the Schrodinger equation and the Helmholtz 
wave equation [131ll33ri35j in order to introduce an index of refraction for matter waves. 

For the case of a constant, complex potential, 

yopt = -2^n2^(/o(fiK))^, (6.10) 
the stationary Schrodinger equation takes the form 

OAJ 

W(r) + -^[S-Kpt]V'(r) = 0. (6.11) 

By analogy with the Helmholtz equation one can introduce a wave number associated to the propagation of 
the scalar field ip (r) with energy E, 
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2M 

^'=y-^[i?-Kpt]. (6.12) 

The effect of the optical potential on the dynamics can thus be equivalently described by an effective index 
of refraction 



- - K - t-^^ (6.13) 



where K = ^rXMEjfi is the wave number associated to the de Broglie wavelength of the tracer particle. 
Since the optical potential can be non-Hermitian in general, the index of refraction is in general complex. 

In the situations described above, where the last, incoherent term of Eq. (|6.9p can be neglected, the matter 
wave dynamics can thus be effectively described by a wave equation with the optical potential Eq. (|6.8p . 
The beam state may then be equivalently expressed as a wave propagation in a medium characterized by 
the refractive index (|6.13p . in complete analogy to the optical case. 

For dilute gases the complex index of refraction is very close to unity, such that 



l + 2-^^(/o(/^K)),. (6.14) 



Its real part 



ni = l + 27r^— Re(/o(;iK))^, (6.15) 



describes the phase shift due to the background gas, while its imaginary part 

n2 = 2^^— Im(/o(ftK))^ (6.16) 



2hK'^lM ' 

determined by the loss term of the classical linear Boltzmann equation, accounts for the attenuation of the 
coherent beam. 

We stress that this analysis in terms of a refractive index relies on neglecting the last term in Eq. (|6.9p . 
This is only valid under very particular experimental conditions. These might be realized in a Mach-Zendcr 
interferometer provided just one of the spatially separated paths interacts with a background gas and 
particles whose momentum has changed due to this interaction arc blocked by the interferometer apertures. 
In this case only the part of the incoming beam scattered in the forward direction contributes to the signal 
collected at the outlet of the interferometer. Note in particular that the preservation of the trace implies 
that the loss and gain term do have the same weight, such that one is not allowed to neglect the last term 
of Eq. (|6.9[) because of its smallness with respect to the anticommutator term. 

At variance with the typical intuition developed for electromagnetic waves the imaginary part of the index 
of refraction, which describes the exponential decay of the beam intensity, does not involve true absorption. 
The tracer particles are obviously not absorbed by the background gas. but the coherent coUimated part of 
the beam entering the interferometer and responsible for the final interference pattern gets reduced. This 
is the reason why the reduction of intensity in the signal does not lead to a corresponding reduction of the 
visibility as in the case of coUisional decoherence [39]. Over a distance L travelled with velocity V = P/M 
the reduction factor is given by the exponential exp (— M°"' (P) L/V). Comparing this quantity with the 
damped intensity cxp {—2n2KL) of an electromagnetic wave with wavenumber K in a medium characterized 
by an index of refraction with imaginary part rt2 gives indeed back Eq. (|6.16p . 

Specializing to a rotationally invariant scattering amplitude, /(p/,Pi) = /(i?;-Erci), with d the angle 
between incoming and outgoing relative momentum, and E^^-j ~ pf/2m^: the kinetic energy in the center of 
mass, one can explicitly evaluate the thermal average introduced in Eq. 
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Here V = P/M is the velocity of the interfering particle and a Maxwell-Boltzmann distribution (|2.6p is 
assumed with Vfj ~ pp/m. 

This expression, derived here from the quantum linear Boltzmann equation [32], agrees with the formula 
in [139|140j . which was recently obtained along a very different line of thought. While it has been used in 
the most recent interferometric experiments with lithium atoms [38], it deviates from previous theoretical 
results |141I142| used for the analysis of the earlier experiments |37I136I1^ . The difference lies in the 
thermal averaging procedure employed, and the averages in Eqs. (|6.ip . (|6.17p . which follow naturally from 
the linear Boltzmann equation, should be considered correct, as also pointed out in [140| . We note that still 
different methods of deriving a gas induced index of refraction were recently presented in |79|144j . 

It was already emphasized that the notion of a gas induced index of refraction can only be applied to a 
restricted set of experimental situations. One may describe phase shift and particle loss phenomena with 
this concept, but it does not apply to proper decoherence phenomena, where the purity of the motional state 
in the beams gets dynamically reduced due to the interaction with a background gas. As described in the 
following section, we have to return to a density matrix description in order to account for such effects. 

7. Collisional decoherence of matter waves 

The concept of environmental decoherence |145|146|147|148|149|29j is often drawn upon to explain how 
the emergence of classical behavior comes about if physical systems get larger and more complex. Though 
formulated in the same framework of open quantum systems also used to describe friction and thermalization, 
decoherence should be distinguished from such dissipative phenomena since it usually takes place on much 
shorter time scales and even in the absence of relaxation processes with a classical analog. 

The main point is to acknowledge the crucial role played by the quantum correlations arising between 
the system and environmental degrees of freedom. The combined state of system and environment is no 
longer separable after the interaction, resulting in a loss of purity in the reduced quantum state, which is 
obtained by disregarding the environmental state by means of a partial trace. In a complementary view, 
which highlights the quantum aspects of the correlation, one may view the interaction as constituting an 
information transfer from the system to the environment. 

In our case of a quantum test particle, collisions with a background gas have the effect of unsharp position 
measurements because the scattered state of a gas particle carries away some information on the position of 
the tracer which would be available, in principle, to an outside observer. One therefore expects that if the 
test particle is initially in a superposition of two positions, as in an interferometer, the coherence of this state 
will be reduced the more the better a typical scattering event can distinguish between the two positions. 
Once the coherence expressed by the off-diagonal elements in the position representation has vanished, the 
motional state will be indistinguishable from a classical mixture. 

In the following we discuss how the dynamics of such collisional decoherence phenomena are described 
by the quantum linear Boltzmann equation. We start with the case of pure collisional decoherence, which is 
present even in the limit of an infinitely massive particle, where no dissipation takes place. We also provide 
the quantitative connection to recent interference experiments [341 , where decoherence due to a background 
gas was observed as a loss of visibility in the quantum interference fringes. After that, we discuss the related 
topics of the loss of coherence in case of momentum superpositions and in the Brownian motion limit. 

7.1. The limit of a massive tracer particle 

The incoherent part of the quantum linear Boltzmann equation (|2.2p simplifies if the test particle is 
much more massive than the gas particles. The function L from Eq. (j2.3[) . which provides the details of the 
collisional interaction, is then replaced by 
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Tl T> /"-gas,/ Q , ^\ f ,Q 

L(p,P;Q) ^ J_2_/|^p^Q__,p^Q + _j ^^,|^p^Q + _ 



(7.1) 



Hence, in this limit the Lindblad operators no longer depend on the test particle momentum P. One can 
now perform the dkj^-integration in Eq. (|2.2|) . such that, recalling Eq. (|2.8p and Eq. (|2.9[) . one gets the much 
simpler expression 



dQAff„i(Q;Q) 



p_c'QX/n^^-,Q.x/n 



(7.2) 



Here the classical rate M^^ (Q; Q) must be evaluated with m» — ^ m. It gives the rate for a particle with 
vanishing momentum P = to experience a momentum change of Q. 

It is now helpful to introduce the probability distribution of the possible momentum transfers according 

to 



P(Q): 



M^l (Q; Q) 



(7.3) 



/dQ'A/^|(Q';Q')' 

The normalization, i.e., the total rate of collisions is given by the classical loss rate known from Eq. (|2.12p 



:Af„t (0) = / dQ'M^,'(Q';Q'). 



(7.4) 



It is determined by the total scattering cross-section and therefore always finite. 
We can thus write Eq. (|7.2|) in the compact form [15] 



dQr{Q)(p 



It puts into evidence that the incoherent dynamics is simply given by momentum kicks 
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(7.5) 



(7.6) 



characterized by a random momentum transfer Q. which is distributed according to Eq. (|7.3p . As discussed 
below, Eq. (|7.5p induces an exponential decay of the off-diagonal elements in position representation, i.e., a 
"localization" of a spatial superposition into a statistical mixture. 

A useful explicit expression is obtained from Eq. (|7.2p by introducing the variables Pi = pUi — p+Q/2 
and p/ = prif = p-Q/2, 



p — e«(Pi-p/)-x/?tpg— i(Pi-p/)-x/?i 



(7.7) 



where P (p) = Airp^ pLjs (pn) is the Maxwell distribution of the momentum magnitudes in the gas. 

Equation (|7.7p corresponds to the result of Gallis and Fleming [13jj once amended for an incorrect factor 
of 27r. This little flaw, which is however experimentally relevant, was not due to a simple calculational 
error, but goes back to an incorrect treatment of a "square of (5- functions" in |40j . Such a strictly ill-defined 
expression arises if asymptotic scattering theory is naively incorporated into a dynamic description. One 
can fix the factor either by circumventing the S-matrix in a weak-coupling treatment [46„150„54.49i77j or by 
formally modifying the S-matrix such that asymptotically outgoing states are not transformed [54166132] . 

The equation of pure coUisional dccohcrencc discussed so far was obtained by the naive limit M ^ oo 
of the quantum linear Boltzmann equation. As a result, the incoherent part Eq. (|2.2p depends only on the 
position operator X of the test particle. Also in the coherent part of Eq. (|2.ip the dependence on P gets lost 
in this limit, since the kinetic energy P^/ (2M) vanishes, while the gas induced energy shift (P) turns into 
an (unobservable) constant. It follows that the dynamics of the momentum distribution of the test particle 
state remains frozen under Eq. (|7.7p . 
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The full quantum linear Boltzmann equation provides of course a much more refined description of the 
dynamics. However, for a massive test particle and at short time scales, where the effects of an energy 
exchange between particle and gas are negligible, an improved version of Eq. (j7.2[) may already be sufficient 
in many situations. A natural improvement is suggested by the fact that the loss term (|7.4|) is evaluated for 
a test particle at rest. 

Provided that the test particle is very massive, such that its momentum varies slowly on the short time scale 
of decoherence, it may be allowed to replace the momentum operator in Eq. (|2.2p by its initial expectation 
value Po = (P)po- This leads to an equation very similar to Eq. (|7.2p . 



Cp^- / dQ M^l (Po + Q; Q) \p - e''J-^/'^pe-''i-^/''l . (7.8) 



Introducing the normalized probability density 

A//;;(Po + Q;Q) 
/dQ'A/>;|(Po + Q';Q') 



<(Po + Q;Q) 



one arrives at 

Cp=-Ttot (Po) j dQPp„ (Q) [p- e^'^-^^/V''-^/''] , (7.10) 

where the rate Ftot (Po) = ^^'^oit (Po) gives the classical rate of collisions experienced by a test particle with 
momentum Po. 

The difference with respect to Eq. (|7.5p lies in the residual, parametric dependence on the momentum of 
the test particle Pq, leading to a more realistic estimate for the total scattering rate Ftot (Po) and thus of 
collisional decoherence. This is of relevance for the quantitative description of experiments on gas-induced 
decoherence |39|55|151|106] . 

It is clear that Eq. (|7.10p cannot account for the long time dynamics of momentum and kinetic energy, since 
the momentum of the test particle here only appears as a classical label. As far as the short time dynamics is 
concerned, we note that the corresponding random momentum transfers ()7.6p are here incorporated in full, 
without the assumption \Ql^X/h\ <^ 1, used in Sect. [5] when considering the limit of quantum Brownian 
motion. 

7.2. Decoherence in position 

7.2.1. The exact solution of pure decoherence 

In position representation Eq. (j7.5p takes the simple form 

{X\Cp\X') = -Ftot [1 - $P (X - X')] (X|p|X'), (7.11) 

where ^-p (X) denotes the Fourier transform of the probability density ()7.3p of momentum transfers, i.e., its 
characteristic function jl52j 

<i>-p{S) = J dQP{Q)c'^-^/^. (7.12) 

The solution of the master equation dp/dt = Cp then shows an exponential decay of the position off-diagonal 
elements, 

(X|pt|X') =e-^'°'[i-*^(^-^')]'(X|po|X'). (7.13) 

The position density is not affected, (X|pt|X) = (X|po|X), since $73 (0) = 1. 

The function $73 (X — X') describes the effect of a single collision event on the statistical operator in 
the position representation, (X|/9|X') $p (X — X') (X|p|X'). It is sometimes called the "decoherence 
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function" since it amounts to the overlap of two scattered gas particle states displaced by X 
quantifies to what extent the gas can distinguish between the two positions. 

It follows from Eq. (j7.3p that in terms of the scattering amplitude it is explicitly given by 



(S) : 
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(7.14) 



In the case of isotropic scattering, the decoherence effect depends only on the distance S = |X — X'| of the 
positions considered. It then takes the form 



Ttot J m J_i 



(cos 61) 



smc I 2p sin ( - j ^ 



(7.15) 



with 9 the scattering angle and sine (x) — sin (x) jx. Due to the infinite test mass the momentum transfer 
is given by the expression Q = 2p sin (0/2), which appears in the argument of the sinc-function. The effect 
of this function in Eq. (|7.15p is thus to suppress position coherences the more the better the scattered state 
with its resolving capacity hjQ can distinguish between two possible scattering positions with distance 5*. 

The second ingredient to Eq. (|7.13p is the total rate of collisions Ftot defined in Eq. (|7.4|) . The effect and 
the rate of collisions are thus fixed by <I>-p and Ftot separately, which can be of great importance in physical 
applications. We note in particular that the function 



(7.16) 



appearing in Eq. (|7.13p is itself a characteristic function, namely that of a compound Poisson process 
|153I152| . Such a process differs from a simple Poisson process in that the jump events, here the collisions, 
do not have a fixed size, while the waiting time distribution between the events is Poissonian. Rather, the 
jump size is itself a random variable, here the momentum transfer in a single collision, which is distributed 
according to V (Q). 

This physical picture behind the coUisional dynamics is best seen by expanding the exponential in 
Eq. (|7.13p . which corresponds to a jump expansion of the solution of the master equation 



oo 

(X|pi|X') = P" W (X - X') (X|po|X'). 



(7.17) 



Here Pn{t) = e~^'°'* (PtotO" /"■■ is the probability to have n collisions during the time t, given by the 
Poisson distribution with parameter Ptot- The associated effect is the multiplication of the initial statistical 
operator in the position representation with the 7i-th power of the decoherence function $p for the effect of 
a single collision. 

Equation (|7.13p is a particular case of a more general class of solutions of master equations for the 
description of decoherence, applying in the presence of translational invariance. In the general case, the 
characteristic function ^' (X, t) corresponding to Eq. (|7.16p is the characteristic function of a Levy pro- 
cess, of which compound Poisson processes are a particular example [29 152] . This result can be obtained 
considering the general structure of translation-covariant Markovian master equations obtained by Holevo, 
and considering its limiting expression relevant for the description of decoherence. This provides a common 
theoretical framework for the description of quite different decoherence experiments | 154 | . The solution of 
such master equations reads 



(X|pt|X')=*(X-X',t) (X|po|X') 



(7.18) 



and due to the general properties of characteristic functions it naturally describes decoherence effects. In 
fact, it is a defining property of characteristic functions that ^' (0,i) = 1, while I'i' (X — X',t) | ^ 1, such 
that diagonal matrix elements are not affected, while off-diagonal matrix elements are generally suppressed. 
In particular, for the characteristic function of a Levy process one has limt^oo 5* (X — X', i) = for X 7^ X' 
|153j . such that all coherences are completely lost with elapsing time. For fixed interaction time t, on the 
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other hand, a variety of position dependencies are possible, from an exponential loss of coherence to revivals 
ending up with a constant decoherence rate. All such behaviors are comprised in the general expression of 
the characteristic function of a Levy process, provided by the Lcvy-Khintchine formula, one of the beautiful 
results of probability theory [152j . 

7.2.2. Experimental tests 

One of the most important experimental verifications of the quantum linear Boltzmann equation comes 
from its application to the quantitative explanation of coUisional decoherence observed in near-field inter- 
ferometry experiments with large molecules |39|55| . In such experiments a velocity selected molecular beam 
from a thermal source passes an interferometer made up of material gratings. After flooding the vacuum 
chamber with various background gases one records the reduction of visibility of the interference fringes as 
a function of the gas pressure. This decrease of the interference contrast is directly related to the loss of 
coherence in the position off-diagonal elements in the state of motion, which was discussed in the previous 
section. 

An important requirement for the experimental observation of coUisional decoherence is that the interfering 
molecule is not removed from the interferometer by the interaction with the background gas. As described 
in Sect. [6l such a loss process would lead to a reduction of the count rate, but it would not diminish 
the interference contrast of the detected molecules. This is the reason why relatively heavy C70 fuUerene 
molecules were required to observe the effect in |39|55j . while the Mach-Zehnder-type atom interferometers 
realized so far cannot observe coUisional decoherence with room temperature gases. We note that the loss 
of visibility in a sodium atom interferometer reported in [155] is therefore due to a different mechanism, 
namely the gas- induced noise at the preparation, i.e., the slit coUimation and at the detection stage of the 
interferometer beam. The atomic motion between preparation and detection is there coherent, in contrast, 
such that this gas-induced effect is not described by the theory of Sect. 17.2.11 

Due to the large mass of the interfering particles and the short time of flight in [39|55j , one is allowed to 
neglect relaxation effects in the quantum description of the interfering test particle. The limit of Sect. 17.11 
thus applies, and one expects in particular that the dynamics is correctly described by Eq. (j7.10p . which 
accounts for the Unite velocity of the molecular beam with respect to the thermal gas. In the position 
representation the equation reads as 



where Pq is now the mean momentum of the incoming fuUerenes. It determines the rate Ftot (Po) = 
A/q^^. (Pq), as well as the decoherence function 



The classical rate Af^J, deUncd in Eq. (|2.10p . includes the differential scattering cross-section for collisions 
between fuUerene and gas molecules. 

The interaction between the gas molecules and the fuUerenes is well described in the experiment by the 
isotropic part of the London dispersion force. It corresponds to a potential of the form U (r) = —C^/r^, 
with an interaction strength Cg, which can be calculated for different gases |156|39j . The associated total 
scattering cross-section can be obtained in a semiclassical calculation [157] ; it depends on the relative velocity 
as 



A quantitative description of the general loss of coherence in a near fleld interferometer due to master 
equations of the form (|7.10p was presented in [151j . However, the situation is somewhat simpler in the 
experiment described in |39I55| since the background gases are at room temperature. This implies that for 
practically all positions in the interferometer the collisions are characterized by large momentum transfers 



{X\£p\X') = -Ftot (Po) [1 - ^v^„ (X - X')] (X|p|X'), 



(7.19) 





(7.20) 
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Q in the sense that Q 3> ft/S', if S denotes the separation of the interference paths. As foUows from Eq. 
(|7.15[) ; aheady a single eoUision then sufhees to completely reduce the off-diagonal elements of the position 
representation, (f>-p (S" 7^ 0) ~ 0, while <i>-p (0) = 1. Only those molecules then contribute to the contrast of 
the interference pattern which experienced no collisions, while the others contribute only to the unstructured 
background. 

In this situation the loss of visibility is determined only by the collision rate Ftot (Po) in Eq. (|7.19p . 
determined by the loss term Ftot (Po) = -^out (Po) iii tl^c hnear Boltzmann equation |106I79| . according to 
the formula [39] 

V = Voc"^'°'(^«)*. (7.21) 



Here Vo is the visibility in the absence of the gas and t is the time of flight of the fuUerenes crossing the 
interferometer. 

It is therefore sufficient to evaluate the collision frequency Ftot (Po)j which is best described by an effective 
scattering cross-section 

Ttot (Po) = ngasJtTeff (Po) . (7.22) 

The expression tToft [Po] for a total cross-section of the form Eq. (|7.20[) has already been presented in 
Sect. 14.3.31 According to Eq. (|4.47p . evaluated for a = —2/5, we have 

47rF(9/10) (i^C^\^" I ( 3 3 f Po V] 

rtot(Po) = n,as^-^(^^j -p^F^^-Wr-yjhTj )■ ^^-2^) 

In the experiment, the beam velocity is small compared to the thermal velocities of the gas, Pq/ (Mvp) <^ 1. 
The expression can then be approximated as ^39.55. 106j 



47rF(9/10) /37rC, 
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(7.24) 



One can introduce a critical gas pressure for the onset of decoherence in terms of the effective scattering 
cross-section introduced in Eq. (|7.22p , 

MkBT 

Po = p 7-25) 

One therefore expects an exponential loss of visibility of the interference fringes with growing gas pressure 
according to 

V = V^e-P^P°. (7.26) 

This behavior was observed in the experiments |39I55| , in quantitative agreement with theory, see Fig. 17.11 

Finally, we mention that the reduction of interference contrast due to momentum exchange effects was 
the subject of several other experimental p?8.159 160 161.162.163] and theoretical |151I164 154.165il66.167] 
investigations. While the physical mechanism is here the scattering, absorption, or emission of photons, the 
dynamic effect on the motional state of the interfering particle can still be described with a master equation 
of the form (|7.5p . (|7.1ip . All that needs to be done is to adjust the rate of events Ftot according to the 
physical mechanism considered, as well as the decoherence function $73 for the effect of a single event |154] . 
Studies of other environmental decoherence effects in interferometers were presented in |168I169I170I171| . 
These are based on master equations characteristic of linear coupling models, which are of the form of the 
Brownian motion master equation discussed in Sect. 17.41 
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pressure (in 10"^ mbar) 

Fig. 7.1. Experimentally observed visibility of the interference fringes of C70 fuUerenes in an argon gas at room temperature, as 
a function of the gas pressure (dots) . The solid line gives the prediction of decoherence theory. It is based on Eq. 117.261 1 , but also 
accounts for gas-induced modifications of the particular gravitational velocity selection scheme employed in the experiment. 
(Figure taken from 1551 ) 

7.3. Decoherence in momentum 

As described in the previous section, it is relatively simple to evaluate and to experimentally test the decay 
of coherence in superpositions of position eigenstates. The analogous case of superpositions of momentum 
eigenstates is more involved. On the theoretical side, the momentum operator of the test particle can no 
longer be considered a classical label in the limiting expression of the quantum linear Boltzmann equation, 
but one deals with a master equation for the operators X and P of the test particle. From the experimental 
point of view, it is difficult to conceive experiments testing decoherence in momentum space. The source 
must generate a state showing longitudinal momentum coherences, which requires a non-stationary beam 
preparation, as well as an interferometric setup permitting to monitor the presence of such coherences 
[1721173] . This is complicated by the fact that the decoherence in position discussed in Sect. 17.21 is always 
present and may quickly dominate. It affects superposition states of different momenta, because the free 
evolution of such states unavoidably leads to a superposition of spatially separated wave packets. 

In order to study the decoherence dynamics of momentum superpositions we exploit the quantum trajec- 
tory algorithm presented in Sect. 14.3.21 taking as initial condition a coherent superposition of momentum 
eigenstates. We start with an initial state of the form 

N 

|^(0))=^a,(0)|UaO)), (7.27) 

where the complex amplitudes Ui satisfy the normalization condition X^i^i l^^i (O)l^ — 1- discussed in 
Sect. 14.3.21 the state at a later time t will still be in a superposition of N momentum eigenstates, but with 
different amplitudes ai (t) and scaled momenta (t) 

N 

|^(i))=5]a.(t)|U,(i)). (7.28) 

i=l 

Given a jump into the state Eq. (|7.28p occurring at time t, the deterministic time evolution described by 
Eq. (|4.39p only modifies the amplitudes. 
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^-iE{U^{t))T/h^-^{Ui{t))r/2 

a^ {t + T)= a, (i) , (7.29) 



with E (U) ~ MvpU"^ /2 the kmctic energy. The cumulative distribution function F (r) of Eq. (|4.40p for the 
waiting time until the next jump is now given by a sum of exponential functions, 

N 

i^(r) = l-^K-(i)|'e-f(^^(*»^ (7.30) 



U,(i)^U,(t) + ^K (7.31) 



The jump modifies amplitudes and momenta as 

U,(i)-U,(t) + ^ 
cti (t + r)^ f^a^ {t + t) 

with 

-(K/2+K-Ui(t))V2 

/.= , . (7.32) 

v/Ef=il«.Wl'e-(WK.u,(t))^ 

While all momentum eigenstates in the superposition are shifted by the same amount, their different am- 
plitudes are modified in different ways. The momentum transfer K characterizing the jump is distributed 
according to Eq. ()4.42p with the state I?/" (i)) given by Eq. (|7.28p . The probability density Eq. (|4.55p is 
therefore replaced by the mixture 

N 

^(A',e)=^A,7'(A-,e^), (7.33) 

i=l 

where the weights A,; are given by 

^ |a,(t)|^e-r(^-W)"f(^,(t)) 
' Ef=il«.Wl'e-f(^^(*))-f(C/,W)- 

The probability density V {K,£) is defined in Eq. (|4.55p and = • K/ (UiK) is the cosine of the angle 
between K and XJi (t). 

We now specialize to a superposition of two momentum eigenstates with equal amplitudes ai (0) = 
a2 (0) = 1/V2, 

|V'(0)) = -i=(|Ui(0)) + |U2(0))), (7.34) 

and evaluate the temporal behavior of the matrix elements in position representation. The off-diagonal 
elements provide a measure for the coherence in the state of motion 



C (t) = E 



"1 (0"2 (0 



ai (0) a* (0) 



(7.35) 



This quantity is plotted in Fig. 17.21 for the choice Ui (0) = — U2 (0) = Uq, such that the distance in 
momentum space is given by |Ui (0) — U2 (0)| = ^Uq. One observes an exponential decay in time in analogy 
to the behavior (|7.13p of the position coherences. 

Based on this observation we now assume that the decay has the generic form 

C(t) = e-^(^")*. (7.36) 

We proceed to evaluate how the decay constant A depends on the value Uq characterizing the initial mo- 
mentum separation and to compare the analytic result with the simulations. 
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0.2 0.4 0.6 0.8 1 1.2 

Fig. 7.2. Semilogarithmic plot of the coherence C (t) of a superposition of momentum eigenstates, as defined in Eqs. 17.3411 . 
1 17. 351 1. One observes a clear exponential decay, in analogy to the behavior of position superpositions. The plot corresponds to 
a constant cross-section with Uq = A and a mass ratio of m/M = 1; the time is given in units of the inverse total scattering 
rate g^Hll. 

The quantum trajectory method presented in Sect. 14.3.21 is equivalent to solving the master equation. 
One can thus assess the decay constant by noting that for short times C (t) either does not change, with a 
probability 1 — r{Uo)t, or it becomes /1/2 if a jump takes place. Neglecting higher order jump processes 
one thus has 

C (t) « 1 - A (f/o) t « 1 - f (C/o) t + f (C/o) t(/i/2), (7.37) 

where (■ ■ • ) denotes the average over the different possible momentum transfers. Exploiting the expression 
Eq. (|7.32p and averaging with respect to Eq. (|7.33p one comes to 

A(C/o)«f(C/o)(l-(/i/2)). (7.38) 
For the case of a constant scattering cross-section, it takes the explicit form 

A(C/o) = f(C/„)~r0^^^^, (7.39) 

with defined in Eq. (|4.28p . As shown in Fig. 17.31 this prediction is in very good agreement with the values 
of A(Uo) extracted from the numerical simulations for various initial values of the momentum Uq. Similar 
results can be obtained for different scattering cross-sections, the degree of accuracy of the analytic formula 
depending on the considered scattering cross-section. 

A direct comparison of relaxation and decoherence rate in momentum space can be made for the case 
m/M <^ 1 and assuming that the momentum separation is large compared to the thermal momentum of the 
gas L/q ^ 1. The relaxation rate is then given by ry as in Eq. (|4.33p . while the decoherence rate Eq. (|7.39p 
can be approximated for Uq ^ I hy A (Uq) w TpUo. The ratio of relaxation and decoherence rate 

rj 8 m 

indicates that momentum decoherence takes place on a much shorter time scale than relaxation effects. 

The study of how decoherence takes place in momentum space is relevant also for the description of the 
transition from the quantum linear Boltzmann equation to the classical linear Boltzmann equation, already 
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Fig. 7.3. Plot of the decoherence rate A in units of the collision rate as a function of the initial scaled momentum Uq. The 
points represent least square fits of the simulated data for the mass m/M = 1, while the line represents the analytic estimate 
given by Eq. ll09t . 

discussed at the end of Sect. 14.21 As observed in Sect. 12.21 the diagonal matrix elements of the quantum 
linear Boltzmann equation coincide with the expression of the classical linear Boltzmanu equation, provided 
one considers the quantum scattering cross-section. The dynamical suppression of the off-diagonal matrix 
elements in the momentum representation as a consequence of decoherence thus amounts to the onset of 
classical dynamics as far as the momentum observable is concerned. 

7.4. Decoherence in the quantum Brownian motion limit 

Having discussed how the full quantum linear Boltzmann equation can be used to describe decoher- 
ence both in position and momentum space, let us now address its description in the Brownian motion 
limit of the equation. The resulting master equation, which corresponds to approximating the environ- 
ment by a bath of linearly coupled harmonic oscillators, is most often used to discuss decoherence effects 
|145ll46ll47ll48ll49l29j . In the simplest case, it leads to an exponential decay rate of the off-diagonal ma- 
trix elements in the position representation, which depends like a Gaussian on the position separation. 
This picture corresponds to the common lore, which however obviously cannot cover all interesting physical 
situations, neither theoretically |174I175I176TT77] . nor experimentally [158I159I160I161I39I162] . 

We start with the master equation obtained by neglecting the friction term in the equation (|5.15p of 
quantum Brownian motion. 

>Cp = - ^ [Ho, P] - ^ E [X^ ' [X^"^]] - ^ E ' [P^' • (7.41) 

Since the friction term in Eq. (|5.15p accounts for energy dissipation, one expects that Eq. (|7.4ip still describes 
the short time decoherence phenomena associated with quantum Brownian motion. 

7.4.1. Position representation 

In the position representation, the general solution of Eq. (|7.4ip can be most easily expressed in terms of 
the solution of the free Schrodinger equation |118ll78ll79j . Denoting the freely evolved state at time t as 
pf , we have 
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X') t 



[Dpp/AM^) e 



\^ XX 

Z/2 



iDpp/3RP)t^] 



Y • (X - X') t 



h2M [D,, + iDpp/3R'P)t^] 



X exp 



y2 



(X-Y|pf|X'-Y) 



(7.42) 



It is easy to recover from Eq. (|7.42p the standard result of elementary decohercncc models, according to 
which the loss of coherence in position space is governed by a Gaussian function in the distance |X — X'| 
[180I40| . In the limit of a very massive particle and short times, the spreading due to the Hamiltonian term 
can be neglected, [Dpp/3M^) ~ 0. One is then left with a convolution of the initial statistical operator in 
position representation with a Gaussian of width 2Dxxt- On short time scales, and for high temperatures, 
the Gaussian is strongly peaked. One is then left with the expected expression 



(X|p,|X')=exp 



D 



^(X- 



X')'0 (X|po|X') 



(7.43) 



7.4.2. Momentum representation 

In the momentum representation, it is more convenient to directly express the solution with reference to 
the initial state po 



(P|p*|P')=exp 



S2 



/ 1 



(P~P') t exp 



3/2 



\ 4:7rDppt 



dQ exp 



_]_Dpj 
12 ft2 



P-P' 



M 



t ] t \ exp 



h 



P'' 



2M 



ADppt 



P-Q|po|P'-Q) 



(7.44) 



Looking at the position matrix elements Eq. (|7.42|) . one notes that the quantum position diffusion term with 
coefficient Dxx mainly corrects the detailed behavior in the exponential suppression of position coherences. 
In the momentum representation, one has two distinct decoherence mechanisms instead. 

The first factor at the right-hand side of Eq. (fr44|) is just due to the i^j^ i^j^P]] term in Eq. (fTiTj) . 

in perfect analogy with the X)j=i [^jj[^jiP]] term causing decoherence in position. The remaining part 
comes from the joint effect of free evolution and the term causing decoherence in position. As already 
mentioned, initial superpositions of different momenta lead to position superpositions as time evolves, such 
that localization effects also destroy coherences in momentum representation. This asymmetry between 
position and momentum is clearly due to the presence of the kinetic term. Indeed, Eq. (|7.42p and Eq. (|7.44p 
acquire the same functional form in the limit A/ — > oo, which removes the Hamiltonian term from the master 
equation. 



7.4.3. Relevance of the quantum correction term 

We now estimate the relative weight of the two contributions at the right-hand side of Eq. (|7.44p . so as 
to gain an insight on the relevance of the quantum correction term with coefficient D^x ■ In order to do that 
we recall the expressions (|5.18p and (j5.19p for the diffusion coefficients in terms of the characteristic thermal 
scales. One thus obtains an expression involving the friction constant 77, as well as the temperature and the 
mass of the test particle. 



(P|pt|P')=exp 



/3 (P - P') 
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\4:Tr Mrjt 
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dQ exp 
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Q|po|P'-Q). 



(7.45) 
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This equation tells us that momentum coherences arc exponentially suppressed for P — P' exceeding the 
thermal momentum spread {P'^)/3 = 3M//3. The first factor corresponding to the quantum term, however, 
is the dominating contribution only for times t < f3h, while for longer times the effect due to the spatial 
dispersion dominates and grows very quickly with elapsing time. 

8. The dynamic structure factor 

We now assume a slightly different point of view on the classical and the quantum linear Boltzmann 
equation, by rewriting it in a way which puts into evidence the statistical properties of the gas, as well 
as the role of the energy and momentum transfer. This is achieved by expressing the collision kernel in 
terms of the dynamic structure factor of the background gas, a two-point correlation function of its density 
fluctuations. This quantity, which is popular in the analysis of experiments where microscopic probes scatter 
off macroscopic samples, is defined as the Fourier transform of the time dependent density autocorrelation 
function of the gas with respect to both the energy transfer and the momentum transfer of a single collision. 
Its appearance in the quantum linear Boltzmann equation suggests a natural interpretation of the ensuing 
dynamics as driven by the density fluctuations in the medium, in accordance with Einstein's understanding 
of Brownian motion |181| . We emphasize that the results discussed in this section are valid for a Maxwell- 
Boltzmann distribution for the gas particles, while they generally rely on further assumptions for more 
general distributions. 

8.1. An alternative formulation of the linear Boltzmann equation 

8.1.1. The classical expression 

Let us go back to the classical linear Boltzmann equation written in terms of an explicit energy conservation 
as in Eq. (|3.1[) . In Sect. [3] it serves as a starting point for the heuristic motivation of the quantum Boltzmann 
equation, where a decisive step is the decomposition of the momenta of gas and test particle into components 
parallel and perpendicular to the momentum transfer Q. This permits one to express the energy conserving 
delta function in the classical linear Boltzmann equation in terms of the parallel momentum components only. 
In the scattering cross-section, in contrast, only momentum components perpendicular to the momentum 
transfer Q appear, as can be seen in Eq. (|A.8|) . 

We can now observe that the Maxwell-Boltzmann distribution (|2.6[) is invariant under rotations and that 
it factorizes in Cartesian coordinates. This implies that for any fixed direction Q one can write 



denote probability densities in two and one dimensions, respectively. It follows that also the momentum 
distribution can be split into factors depending either on the parallel or the perpendicular components of 
the gas momenta with respect to the momentum transfer Q. 

Exploiting these facts, one can perform the integral over the parallel component of the gas momentum in 
the linear Boltzmann equation Eq. (|A.8p . This can be done without specifying the scattering cross-section 
since the latter depends only on the orthogonal components of the gas momenta. One thus obtains the 
compact expression 




(8.1) 



where 




(8.2) 



(8.3) 
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9r/(P) = ^ / dQ y'^dk^A*(^'^)(k^)a(^rel(k^,PxQ)-|,rel(k^,PxQ) + §) (8.4) 
X [Smb (Q, P - Q) / (P - Q) - Smb (Q, P) / (P)] , 



upon defining 



Smb (Q,P) = y" dk/.^ (k) S 



2M 2m 2M 2m 



.[id] 



Qll 



1 



The second line follows from the identity Eq. (|A.4[) and the normalization is fixed in Eq. (|8.3p . As will be 
discussed in Sect. 18.21 the function 5mb (Q, P) is the dynamic structure factor of a gas of Maxwell-Boltzmann 
particles. 

By further introducing an averaged scattering cross-section, involving the two-dimensional distribution 
(j8.2[) over the perpendicular gas momenta, 

cTav (PxQ, Q) = J^^ dki ^if'^ (ki) a (^rel (k^, P^q) - §, rel (k^, P^q) + , (8.6) 

Equation (|8.4p assumes the compact form 

9r"/(P) = ^ / dQ (P±Q, Q) [5mb (Q, P - Q) / (P - Q) - Smb (Q, P) / (P)] • (8.7) 
mi J 

By exploiting the factorized expression of the gas distribution function Eq. (|8.ip we thus obtained an 
expression of the classical linear Boltzmann equation which involves only an integral over the momentum 
transfer Q. Note that if the scattering cross-section only depends on the momentum transfer, as is the case 
in the Born approximation, the kj^-integration can be done even in the absence of a factorization property 
like Eq. (|8.ip . The averaged scattering cross-section (|8.6p then coincides with the proper cross-section. 

The integral in Eq. (|8.5p can be calculated explicitly, leading to an expression which is best stated in 
terms of the momentum Q and the energy E transferred to the test particle, whose momentum prior to the 
collision is P. We take these quantities as positive if the test particle gains momentum or energy. According 
to Eq. ((^TT|) the energy transfer is E (Q, P) = Q^/ {2M) + Q -P/M, such that it depends only on P||q. The 
function Smb can therefore be expressed equivalently as a function of Q and P||q, or of Q and £'(Q,P), 



.„3,.,p,^.„e,Q,.(Q,p,,^;^l»p(-A2™^). ,s.a, 

This exact expression implies in particular 

(p^Q + £^§ + I^PiIq) =4''^ (p^q) ^mb (Q,P||q) , (8.9) 

where we stress explicitly that the function Smb depends on P only through Pyq. Using Eq. (|8.9p one 
immediately finds that Eq. (|8.4p and Eq. (|3.2p are two equivalent forms of the same equation. 



8.1.2. Generalization to the quantum case 

The particular form (j8.4p of the classical linear Boltzmann equation suggests an analogous way of writing 
the quantmir linear Boltzmann equation |58l42l47ll82ll83j . 
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Q 



-2 S^mb(Q,P) 



/ ( rel (k^, P^q) - |, rcl (k^, P^q) + | 



(8.10) 



This is exactly equivalent to Eq. (|2.2p for a gas of Maxwell-Boltzmann particles, described by Eq. (|2.6p . 
as can easily be seen using Eq. (|8.9p . Note that the square root in Eq. (jS.lOp is well defined since Smb is 
positive, which is a general property of the dynamic structure factor, as discussed in the following section. 



8.2. Properties of the dynamic structure factor 



8.2.1. Expression in terms of the density correlation function 

Let us consider a quantum many-body system, e.g. the gas of free particles considered in the previous 
paragraph, described in second quantization by a field operator ijj (x), satisfying canonical commutation or 
anti-commutation rules. One can then consider the Fourier transformation of the operator density n (x) = 
"0* (x) '0 (x) given by 



Pq = y" dx e-'^ ^'/^n (x) . 

The associated spectral function 
1 1 



2TThN 



dt c^'^'/^Pc^p^it)), 



(8.11) 



.12) 



with N the total number of particles, is known as dynamic structure factor [184|185|186j . though other 
names such as "spectral density function" , "scattering function" and "scattering law" are also used in the 
literature. It can be expressed in general as the Fourier transform with respect to energy and momentum 
transfer of the time dependent density autocorrelation function 



as 



G(x, i) = 1 y" dy (n (y) n (x + y, t)) 

5(Q,S) = ^y"diy" dx e'(^*-Q-"'/'*G(x,t). 



.13) 



.14) 



By evaluating the dynamic structure factor for the case of a system of free distinguishable particles one 
obtains, recalling Eq. (|4.1ip for the energy transfer. 



5(Q,i?(Q,P)) 



dp 



np)5 i?(Q,P) 



(P - Q)' 



2m 



2m 



(8.15) 



Here (rip) denotes the expectation value of fip, i.e. the mean number of particles with momentum p. If we 
consider a Maxwell-Boltzmann distribution, such that 



(hp) ngas {27rh) pp (p) , 
one comes back to the expression introduced in Eq. 



.16) 



(see e.g. [47ll87p . 
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It thus appears that the dynamic structure factor Smb, which can be used to express the scattering 
kernel of the linear Boltzmann equation (jS.lOp . has an important and transparent physical meaning. It 
is the Fourier transform of the correlation function of the density fluctuations in the medium, reflecting 
the molecular nature of the gas, which determines the coUisional interaction between test particle and gas 
particles. It thus embodies in a suitable correlation function the original intuition by Einstein which led to 
the correct understanding of Brownian motion: the dynamics of the test particle is driven by collisions taking 
place due to the discrete nature of the medium, whose density fluctuations provide the coupling mechanism. 



8.2.2. Connection to the laboratory-frame scattering cross-section 

The dynamic structure factor is particularly useful in describing the scattering of microscopic probes 
off macroscopic samples, since it allows one to express the energy dependent cross-section in terms of a 
correlation function of the many-body system [188|185j . It was in fact flrst introduced by van Hove in order 
to describe neutron scattering |189j . 

Using Fermi's golden rule and the Born approximation one can show that the laboratory-frame cross- 
section S for the scattering of a test particle off a many-body system is characterized by the dynamic 
structure factor S {Q,E) of the many-body system |131| 

J2y l\,r2 pi j\ r2 pi 

'V{Q)\^S{Cl,E) = ——\fBm''S{Cl,E). (8.17) 



Here a scattering event is considered in which a particle of mass M changes its momentum from P to 
P' = P -)- Q, absorbing an energy E according to (|4.1ip . and V^(Q) is the Fourier transform of the two-body 
interaction potential 

y(Q) = j dxc'^ ^'/^V (x) , (8.18) 

related through Eq. (|2.15p to the scattering amplitude in Born approximation. 

More generally, provided multiple scattering effects can be neglected, one has jl88) 

SiQ.E), (8.19) 



dQp'dEp' dnp. 

where dS/dfip' describes the single scattering event in the laboratory frame. 



8.2.3. Detailed balance and the stationary solution 

For a strongly interacting many-body quantum system the dynamic structure factor is no longer analyt- 
ically tractable as in Eq. (|8.8p . but it can only be measured in suitable scattering experiments. However, 
due to its very definition the dynamic structure factor acquires many interesting general properties. First of 
all, it is always positive as a consequence of Bochner's theorem |152j . because the density autocorrelation 
function Eq. (|8.13p is a positive definite function. This property is to be expected on physical grounds, since 
the dynamic structure factor is directly proportional to a scattering cross-section according to Eq. (|8.17p . 
The dynamic structure factor further obeys the so-called detailed balance condition 

S{Cl,E)^e-'^'^S{-Cl,-E), (8.20) 

which is the crucial property for ensuring the existence of a stationary solution of both the classical and the 
quantum linear Boltzmann equation. In order to confirm this property one has to consider Eq. (|8.12p and to 
observe that /5q = /5_q , as follows directly from the definition Eq. (|8.1ip . and to exploit the fact that the 
expectation value is obtained with respect to the equilibrium state of the many-body system pEQ oc c'^^ , 
where H is the full Hamiltonian. The property Eq. (|8.20p therefore holds generally, for any many-body 
system in a thermal equilibrium state. 

In our case, where the energy transfer can be expressed by Eq. (|4.1ip . the detailed balance condition can 
also be written as 
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S (Q, P) = e-'^^t^^P^^ (-Q, P + Q) 



(8.21) 



It is now simple to prove that Eq. (|8.2ip ensures the existence of a stationary solution of the canonical 
form peq(P) oc exp (— /3P^/ (2A/)) for both the classical and the quantum linear Boltzmann equation. For 
/ (P) peq(P) the expression 



^MB (Q, P - Q) PEQ (P - Q) - Smb (Q, P) Peq (P) 



.22) 



is odd in Q, such that Eq. (|8.4p implies 9™"peq(P) = 0. Similarly, in the quantum case the expression 
Eq. (|8.10p for an operator p = J^(P), which is only a function of the momentum operator P, simply reads 



CHP) 
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dQ / dkj 

X [Smb (Q, P - Q) ^^(P - Q) - Smb (Q, P) HP)] ■ 



/ rel(ki,PiQ) 



|,rel(k^,P^Q) + | 



(8.23) 



For the case i'(P) !-'eq(P) the condition Eq. (|8.2ip again implies that the integrand is an odd function of 
Q, such that £i^EQ(P) = 0. Note that the explicit expression Eq. (|8.8[) of Smb is not relevant for the proof, 
which relies only on the property Eq. (j8.2ip . valid for a generic medium. 

As expected, also in the quantum case the existence and the form of the stationary solution do not depend 
on the scattering amplitude, even though it appears operator-valued. This strong correspondence with the 
classical case is due to the fact that the covariance under translations implies that the generator C leaves 
the algebra generated by the momentum operator invariant, as discussed in Sect. 



8.2.4. Fluctuation- dissipation relationship 

Let us now discuss the origin of the detailed balance property, beyond its direct calculation based on the 
definition of dynamic structure factor in Eq. (|8.12p . To this end, we have to touch upon the fluctuation- 
dissipation theorem. 

The dynamic structure factor can be directly related to the dynamic susceptibility studied in linear 
response theory. Denoting the imaginary part of the dynamic response function as x" E!), one finds |186j 

S{Q,E)^--^x"iQ.E). (8.24) 

Since the imaginary part of the dynamic susceptibility, which describes the dissipative part of the response, 
is by definition an odd function of energy, this relationship ensures that the dynamic structure factor obeys 
the detailed balance condition Eq. (|8.20p . We recall that, contrary to the usual perspective in linear response 
theory, we are here concerned with the reduced dynamics of the test particle, thus taking the momentum 
and energy transferred to the particle as positive. 

Introducing the real correlation function of the gas medium 

^^{Q,t)^^{{pQit),PQ}), (8.25) 

for the fluctuations of the operator pQ given by Eq. (|8.1ip . the fluctuation-dissipation relationship can be 
formulated in terms of the dynamic structure factor as follows jl84l56) 

0+(Q,i) = -i / dE cxp(i^t^coth(i^E] (l-e'^^)S{Q,E). (8.26) 



fiJ-oo \ J V2 , 

This is a consequence of the definition of the dynamic structure factor Eq. (|8.12p as well as Eq. (|8.25p . When 
expressed in terms of the imaginary part of the dynamic susceptibility according to Eq. (j8.24p this identity 
gives back the usual formulation of the fluctuation-dissipation theorem, once again taking as positive the 
energy when transferred to the test particle. 
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In this spirit, the expression (j8.17p for the laboratory- frame cross-section can be viewed as a formulation 
of the fluctuation-dissipation relationship. The test particle experiences dissipation due to energy and mo- 
mentum transfer processes described by the cross-section on the left-hand side of Eq. (|8.17p . which is related 
on the right-hand side to the equilibrium density fluctuations, characterized by the dynamic structure factor. 



8.3. Extension to different reservoirs 



In the previous sections, we considered the explicit expression of the dynamic structure factor for an ideal 
gas described by the Maxwell-Boltzmann distribution. The latter implies that the gas particles are either 
distinguishable or that the gas temperature is sufhciently high such that effects due to their indistinguisha- 
bility can be neglected. We note that the result (|8.8p is the same in the classical and the quantum calculation 

For the case of an ideal gas, one can calculate the exact dynamic structure factor also by taking into 
account the effect of quantum statistics. Starting from the density autocorrelation function Eq. (j8.13p for a 
gas of identical bosons or fermions it takes the form 



5BF(Q,i?(Q,P)) = 



^gas 



dp 



(np)BF(l± (np)BF)<5 rB(Q,P) 



2m 



P 

2m 



(8.27) 



which should be compared to Eq. ()8.15p . Here the upper sign stands for Bose-Einstein and the lower one for 
Fermi-Dirac statistics, respectively. The corresponding distribution functions 



1 



p)BF 



-le'^fef ip 1 



(8.28) 



are characterized by the fugacity z — exp (/3/i). The integral (|8.27p can be done, leading to [17] 



1 2TTm^ 1 



±1 



is.s(3{27ThfQl-e0^(Q'P^ 



log 



1 =F z exp 



fj (Q^-l-2m£:(Q,P))' 

" 8m 



1 =F z exp 



Srn 



.29) 



In the limit of small fugacity z rigas (27r/i^/3/m)'^^^, corresponding to negligible degeneracy, the latter 
expression goes back to Eq. (|8.15p . As one can also check directly, 5bf {Qi E) obeys the detailed balance 
condition (|8.20p for both statistics. 

Derivations of the quantum linear Boltzmann equation in the weak-coupling limit, in which a thermal 
reservoir of identical quantum particles is considered, have been obtained in [58j for the case of a bosonic 
reservoir, and more generally in [42], leading to 



''gas 



dQ|/i3(Q)|' 



e^Q-x/V^BF (Q, P)pVSbf (Q, P)e-*Q->^/'"' - - {^bf (Q, P) , p} 



.30) 



This result is the quantum counterpart of a classical linear Boltzmann equation, which is modified in order to 
account for statistical corrections. The latter can be obtained upon setting /ibf (p) = (^^p)bf/ ^'^gas (27r?i)^^ 
and replacing p.^ (p) in Eq. (|3.ip by ^bf (p) (1 ± "gas {2T:Kf pbf {p')), which leads to 



^r"./(P) = ^ J dP'dp'dp5 



(El 
\2M 



2m ^ 2M 



P 

2m 



(P' + p' - P - p) 



X(T (rel (p, P) , rel (p', P')) Pbf (p') (1 ± "gas {^irhf mbf (p) / (P') 



-MBF (p) (l ± rigas {2TTh:f PBF (p')) / (P)J • (8-31) 

The statistical corrections in this modified classical linear Boltzmann equation are the analog of the modi- 
fications introduced by Uehling and Uhlenbeck to the original classical equation [9]. They express the fact 
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that quantum correlations in the gas can be relevant even in the absence of any self-interactions, just because 
the occupancy of the different momentum states may enhance or suppress the related scattering rates due 
to statistics. 

We mention that a simplified form of the master equation (jS.lOp was used in [58] to study how the effects 
of quantum statistics in a Boltzmann gas affect the reduction of visibility in the interference fringes due to 
collisional decoherence. 

Finally, it is interesting to note that the comparison between Eq. (|8.30p and Eq. (|8.10p suggests a natural 
heuristic approach to extend the quantum linear Boltzmann equation to real, interacting gases. In place of 
the dynamic structure factor of the ideal gas one just uses the relevant expression of the interacting gas. 
Even if an exact analytic evaluation of the dynamic structure factor is in general impossible for a truly 
interacting many-body system, one may rely either on experimental data from scattering experiments or on 
a phenomenological ansatz, constrained by the requirement of the detailed balance condition. This strategy 
is straightforward for the weak-coupling form (|8.30[) of the master equation. However, complications arise 
if the interaction with the gas is to be treated beyond the Born approximation, as in Eq. (|8.10|) . once the 
momentum distribution in the gas no longer factorizes as in Eq. (|8.ip . 

8.4. Inclusion of internal degrees of freedom 

The quantum linear Boltzmann equation provides a realistic description for the dynamics of a massive 
test particle traveling through a dilute gas of free particles, and its predictions in the quantum framework 
have been tested in recent experiments relevant for atom interferometry |37|39|155|38| . It is however natural 
to consider extensions of the equation, allowing one to cope with a wider experimental scenario. 

In the first instance both test particle and gas particles have been considered as pointlike particles, whose 
only relevant degrees of freedom are those connected to the center of mass. This is obviously not necessarily 
the case and one should consider a proper quantum treatment of both internal and center of mass degrees 
of freedom. 

A hybrid equation, describing the internal degrees of freedom within quantum mechanics, and the center 
of mass degrees of freedom as classical variables, was recently introduced under the name "Bloch-Boltzmann 
equation" |190I191| . It combines the master equation for the gas- induced incoherent dynamics of a iV-level 
system of internal states with the classical linear Boltzmann equation. This Bloch-Boltzmann equation 
thus accounts for the possibility of inelastic collisions connecting the various possible scattering channels, 
which can be of relevance e.g. for spectroscopic experiments [192| . transport phenomena, or more refined 
interferometric experiments testing decoherence. 

A microscopic derivation of the effect of the collisional dynamics on the internal degrees of freedom of a 
very massive test particle at rest was considered in |69I65| . A fully quantum description of both internal and 
center of mass degrees of freedom of a quantum test particle interacting through collisions with a gas has 
been considered in [193j , leading to an equation called "quantum Bloch-Boltzmann equation" . From this 
equation all the others can be obtained as limiting situations in which one of the two kinds of degrees of 
freedom is traced over or treated classically. In this framework, a rich phenomenology of dynamical situations 
can be considered, including non-Markovian effects described by means of the so-called generalized Lindblad 
structure [194119511961197] . 

The "quantum Bloch-Boltzmann equation" can be written in Lindblad form in close analogy to the linear 
Boltzmann equation ([2.ip . However, the statistical operator p now acts on the Hilbert space describing both 
internal and center of mass degrees of freedom. The correction to the free Hamiltonian, which is the analog 
of Eq. ([2.5[) , now reads 




(8.32) 
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where ~ \i) {j\ provides a basis of operators in the space of the internal degrees of freedom, and £ij = hui 
fkVj are the possible transition energies. The multi-channel scattering amplitudes, denoted by (p/, Pi) : 
/(Piii ~^ P/j*)i describe a transition from an m-state with labels Pi,j to an out-state with labels p/,i. 
The incoherent part of the master equation can be written analogous to Eq. (|2.2p . It takes the form 



-l{p,Lt (ki,P;Q,£:)i(kj.,P;Q,f)} 



.33) 



where, apart from the momentum integrations, there is an additional sum over a discrete index £, labeling 
the energy transfer to the internal degrees of freedom in case of inelastic scattering. 
The function L, which describes the effect of the single collisions, is given by 

L (p, P; Q, = f^^l (PXQ, P±q) - ^ + 7W[^^^^ ^el (p^Q, P^q) + I + 



£ij = £ 




F ■ 



Q ®Ey, (8.34) 



which should be compared to Eq. (|2.3p . This implies that the multi-channel scattering amplitudes /y and 
the distribution function ^ of the gas momenta appear again operator- valued in ()8.33p . 

Equations (|8.32p - (|8.34|l represent the natural generalization obtained by combining the two master equa- 
tions for the quantum motion of a point particle and for the internal dynamics of a system at rest, respectively. 
It should be emphasized, however, that a physically stringent derivation of Eqs. (|8.32p - (|8.34p is still miss- 
ing. It is a subject of current research whether this can be achieved by means of the monitoring approach 
described in [65j . which was also instrumental in deriving the quantum linear Boltzmann equation. 



9. Conclusion 



In the present report we discussed a Lindblad master equation for the quantum motion of a test particle 
in an ambient gas, arguing that it is the natural quantum counterpart of the classical linear Boltzmann 
equation. This means that the gas-induced phenomena described by the equation are accounted for not 
just in a phcnomcnological sense, but are incorporated in a non-perturbative fashion by means of the exact 
scattering amplitudes. This microscopically realistic description of the individual scattering processes allows 
one to predict experimentally observable properties such as dissipation and decoherence rates directly from 
the microscopic interaction laws. At the same time, it explains why the structure of the equation is more 
complicated than most master equations used for the description of open quantum systems. 

We presented a fairly straightforward, heuristic motivation for the particular form of the equation, referring 
the reader to the literature for a more stringent derivation |32I41| . It is remarkable that already two basic 
requirements serve to fix the master equation almost completely, namely the translation-covariance of its 
Lindblad structure and the compatibility with the classical linear Boltzmann equation for quantum states 
which are indistinguishable from classical ones. The remaining freedom in this heuristic approach can be 
settled by a plausibility argument, which is corroborated by the fact that the various limits that can be taken 
turn the equation into established master equations. In particular, a generalized Caldeira-Leggett master 
equation of quantum Brownian motion is obtained in the diffusive limit, whose Lindblad structure is as close 
as possible to the classical Kramers equation with microscopically defined friction and diffusion constants. 

Some attention was devoted to analyzing the general structural features and symmetry properties of the 
equation, as well as the predicted approach to an equilibrium state. We saw that many of the relaxation 
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and dissipation effects are shared with those of the classical linear Boltzmann equation, while pure quantum 
phenomena, such as collisional decoherence and the gas induced index of refraction for matter waves, are 
genuine results of the quantum version. A remarkable feature of the equation is that it accounts for the 
interplay of all relevant gas-induced phenomena, from the short-time decoherence dynamics to the approach 
towards equilibrium for asymptotically large times. At the same time, it is clearly impossible to state the 
general solution in an explicit form, even in the simplest case of s-wave scattering with a constant scattering 
length. However, we saw that the Lindblad structure of the master equation admits a stochastic unravelling 
of the master equation in terms of quantum trajectories, which is particularly transparent in the momentum 
basis. 

The validity of the quantum linear Boltzmann equation is expected to cease in situations where the medium 
cannot be considered an ideal, homogeneous Maxwell-Boltzmann gas, or where the scattering interaction is 
inelastic. We discussed an alternative formulation of the master equation in terms of the dynamic structure 
factor of the gas. Apart from providing physical insights, it suggests a natural way of incorporating corrections 
due to quantum degeneracies and self-interactions in the gas, at least in the framework of the weak-coupling 
approximation. As for inelastic scattering, we described a natural multi-channel scattering extension of the 
equation, capable to account for the presence of internal degrees of freedom in the test particle. 

A further important extension of the equation concerns the treatment of an inhomogeneous background 
gas, e.g. due to an external potential. In view of the corresponding classical equation, it seems natural to 
incorporate this by means of an operator- valued gas density rigas (X), rendering the transition rate Afjn a 
function of the position operator of the test particle. However, the non-commutativity of the position and 
momentum operators complicates this task considerably. 

Finally, the limit of a dense and strongly self-interacting gas turns the background medium into a liquid. 
Clearly, the interactions with the test particle are then no longer described by statistically independent 
two-particle collisions, leading to a breakdown of the Markov assumption. In spite of many years of research, 
the field of non-Markovian open quantum dynamics is still in its infancy. However, it will certainly be of 
central relevance to cope with such non-Markovian effects in order to describe a truly quantum microscopic 
dynamics. 

The long path leading to the present version of the quantum linear Boltzmann equation suggests that 
most of these improvements are still a long way ahead, and possibly will require a different mathematical 
framework. Nonetheless, they will be of great importance for a quantitative assessment of quantum transport 
and decoherence phenomena. 
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Appendix A. Classical formulation 

A.l. The classical linear Boltzmann equation 

The classical linear Boltzmann equation can be written in quite different ways, which reflects its compli- 
cated structure and the many variables appearing in it [2] . The most compact and perhaps well-known form 
is given by 

9rV(P)= J dp I dnairel (p, P) , rel (p', P')) Jrci (p, P) 

x[^(p')/(P')-/i(p)/(P)], (A.l) 
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where jrci (p, P) = "-gaslrel (p, P) |/m* denotes the current density in the relative motion. In this equation, 
the values of p' and P' are determined implicitly by the conservation of momentum, p' + P' = p + P. and 
energy. |rel (p', P') | = |rel (p, P)| , as well as by the angles of rotation which connect the relative momenta 
appearing in the differential cross-section a. Moreover, we have only considered the collisional term. 

However, it is much more convenient for our purposes to use a more explicit expression. In particular, is 
permits to clearly demonstrate the connection between the classical and the quantum version of the equation. 

Let us start from the basic expression obtained by turning the original, non-linear Boltzmann equation 
into the linear version. By replacing one of the distribution functions / (p) appearing in the bilinear collision 
term with the equilibrium distribution /i (p) of the gas particles, one gets an equation 



n„ 



acoii^(p) ^ ^ j ^p' j ^ ' j ^ s ( — + ^ - — - ^ ] 6^P' + p' - P - p) 
* ■'^ ^ ml J J J \2M 2m 2M 2m J ^ ^ ^' 

Xd (rel (p, P) , rel (p', P')) [^^ (p') ./ (P') - M (p) / (P)] , (A.2) 

which can be brought into the form Eq. (jA.ip . Here the dynamics determined by collisions between test 
particle and gas particles is put into evidence. The collisions are characterized by a differential scattering 
cross-section a (p/,Pi), together with the constraints of energy and momentum conservation expressed by 
the (5-functions. 

Exploiting the momentum conservation and introducing the momentum transfer Q = P' — P, which 
corresponds to the momentum gained by the test particle in a collision, we obtain 



3rV(P) = ^ /dq/ dp 6 



(p + qr , {p-Qf 



2M 2m 2M 2m I 

xa (rel (p, P) , rel (p - Q, P + Q)) (p - Q) / (P + Q) - M (p) / (P)] ■ (A.3) 

It is now convenient to express the ^-function of energy conservation in terms of the momentum transfer 
and of the momentum components parallel and perpendicular to it. Introducing P||q = (P • Q) Q/Q^ and 
P^Q = P — P||Q and further exploiting the definition of relative momenta given in Eq. (j2.4p one has 

/ (P + Q)^ , (p - Q)2 P2 p2 \ / Q2 1 




2M ' 2m 2M 2ml ^ V 2m. " ' ^P^^' ^"^^ ' ' ^^''^ 

which implies 

rel(P||Q,P||Q) = §. (A.5) 
Since the reduced mass m, obeys m^/m + m^:/M ~ 1, one has the useful relation 

rel(p,P) + Q = rel(p + Q,P- Q), (A.6) 
which, together with Eq. (|A.5|) . leads to 

a (rel (p, P) , rel (p, P) - Q) = a (rel (p||q, P||q) + rcl (p^q, P^q) - Q, 

rel (p||Q, P||q) + rel (p^q, P_lq)) 

= a (^rel (p^q, P^q) - ^Q, rel (p^q, P^q) + ^Q^ . (A.7) 



Eq. (jA.Sp now becomes 



?rV(P) = ^ jdCtj dpa(^rel(px,Pxq)-^Q,rel(px,P±q) + ^Q 



y.5[-^-—Q-vc\ (piiq, P||q) ) [m (p - Q) / (P + Q) - M (P) / (P)] ■ (A. 
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By means of the shift of coordinates 

m Q rn 

m, 2 AI " 



(A.9) 



which leaves p_L invariant, the argument of the delta function simplifies to Q-p, such that it can be evaluated 
explicitly leading to 

a.o«^,p,^!tep I ^ I dk^a(^rel(k^.P^a)-3 „l(k^,p^j,) + 3j (A.IO) 

("^ -ii . s (p,„ - Q)) / (p - Q) - . (k.+i;i§ . ^p.,) / ,p, 

This is the explicit expression allowing one to straightforwardly see the connection between the classical and 
the quantum linear Boltzmann equation. Another useful form is given by Eq. (|A.14[) . 

A. 2. Approach to equilibrium 

We now prove that the relative entropy of a solution of the classical linear Boltzmann equation with 
respect to its stationary solution is a monotonically decreasing function, equal to zero if and only if the 
solution is at equilibrium. Let us first consider the explicit derivative of the relative entropy (|4.16p . 



dg(/|/EQ) 

dt 



= JdP ar"/(P) 



1 + log ■ 



/eq (P) 



(A.ll) 



Exploiting the relations Eq. (|2.12[) and Eq. (|2.1ip for the loss and gain rates, the classical linear Boltzmann 
equation (j3.3p reads 



9™"/(P) = J dP' [M"^ (P' ^ P) / (P') - M""^ (P ^ P') / (P)] . (A.12) 

The fact that the transition rates obey the detailed balance condition Eq. (|4.9p implies that the quantity 

W (P, P') = A-r' (P' ^ P) c'^s^ (A.13) 
is symmetric in its arguments. Wc thus obtain an expression 

c}rV(P) ^ j dV'W (P, P') [/eq (P) / (P') - /eq (P') / (P)] , (A.14) 

that puts into evidence that the integrand is a product of functions which are symmetric and antisymmetric 
under the exchange of P and P'. Using this form, the time derivative of the relative entropy therefore reads 
as 



di?(/|/EQ)_l 



= 1 1 dP 1 dP' (P' ^ P) [/eq (P) / (P') - /eq (P') / (P)] 

X [log (/eq (P') / (P)) - log (/eq (P) / (P'))] • (A.15) 



By noting that 

(X - y) (logy- log X)^0, (A. 16) 

where the equal sign holds if and only li X — Y , this proves that the time derivative of the relative entropy 
is equal to zero if and only if one considers the equilibrium solution. 

Moreover, since a™"/(P) = imphes dH (/|/eq) /dt = 0, it follows that dH (/|/eq) /dt = 0, or equiva- 
lently / = /eq, is a necessary condition for a stationary solution. But we already know that / — /eq is a 
sufficient condition for stationarity, such that the stationary solution is unique. 
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A. 3. The friction coefficient of classical Brownian motion 



We now explicitly derive the microscopic expression for the friction coefficient appearing in the classical 
Kramers equation (|5.1|) . Let us start from the classical linear Boltzmann equation written as in Eq. (jA.ip . 
with the constraint |rel (p, P) | = |rel (p', P')| keeping track of energy conservation in each collision. We put 
/ (P) = j/EQ (P) X where i/eq (P) has the form Eq. (|4.7p of a Maxwell-Boltzmann distribution of the 
test particle. Exploiting energy conservation Eq. (jA.ip now takes the more explicit form 



dtf(P) 



^i^eq(P) / dppV/3(p) / d$ / dOsine 
Jo Jo Jo 

X |rel (p, P) \a |rel (p, F)\)[x (P') - x (P)] , 
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2tt 



dip 



d sin i9 



(A.17) 



where ip) denote the polar angles of rel (p', P') with respect to rcl (p, P), and (8, $) the polar angles of 
p with respect to P. 

As discussed in Sect. [5l in the Brownian motion limit one considers a test particle much heavier than 
the gas particles, such that m/M ^ 1, and the test particle is assumed to be close to thermal equilibrium. 
Under these assumptions the momentum transfer in a single collision, given by 



Q = P' - P = - [rcl (p', P') - rel (p, P)] , 



(A.18) 



can be considered small. We therefore expand the function x (P') iii a Taylor series, thus obtaining up to 
second order 



dtf{V) 



i^eq(P) / dppV/3(p) / d$ / dOsine 



dp 



X |rel(p,P)|f7(i9;|rel(p,P)|) 

X E (p' - tIt^ (p) + ^ E (p' - p). (P' - p)^ 



1=1 



dPjdP, 



d i9 sin 1? 



X(P) 



(A.19) 



Thanks to the identity Eq. (jA.lSp and to the fact that the scattering cross-section only depends on the 
angles i9 and 0, one can perform the integration over the azimuthal angle p appearing only in P' — P. It is 
done by exploiting the integrals 



j'\p{K'-K)^{K'~K)^^2n 



sin^ -ffK^S, 



(3 cos 7?- 1)(1 - COS'S) KjK, 



(A.20) 
(A.21) 



which are valid whenever (i?, ip) denote the polar angles of K' with respect to K. 

To obtain Eq. (|A.20[) let us denote by ei, 62, and 63 the unit basis vectors for the Cartesian coordinates, 
such that Ki = e.; • K, and consider another set of basis vectors e'^, e2ande3 such that 63 is in the direction 
of K. One then has 



(K' - K),^ = e, • (K' - K) (A.22) 

= ei ■ [K' sin cos pe'i + K' sin sin (^e'j + {K' cos -d — K) e';^] . 

Equation (|A.20p is obtained by integrating the azimuthal angle p and recalling K = Kg'^- A similar, but 
much longer calculation leads to Eq. (jA.21[) . 

We are left with contributions in the integrand proportional to rel (p, P)^ and rel (p, P)^ rel (p, P)^, such 
that the integrals (|A.20[) and (|A.2ip can be applied once more to perform the integrations over <i>. The result 
of the integration over the azimuthal angles reads 
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9J(P) = 47r2^iyEQ(P) / dpp^i^inip) dOsine / d?9siii?9 (1 - cosi?) |rel (p, P) 

Jq Jo Jo 

' M p 



xa(,?;|rel(p,P)|) 



■ cos 



i—1 



(A.23) 
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To proceed we exploit the fact that the test particle is assumed to be close to equilibrium, such that 
p/P « m/M , and m/M ^ 1. We expand the expression of |rel (p, P) | in this small ratio, thus obtaining 
up to order yJm/M 



|rel(p,P)|«p 



1 ™ ^ r\ 
1 — — — cos H 



M p 

and therefore in particular 

|rel(p,P)|a(79;|rel(p,P)|)«p 



a{d-p) 



Tft P f d 

— — cos e [a p) + p—a {d]p) 



(A.24) 



(A.25) 



As a result of this approximation one can perform the integral over O, without explicit knowledge of the 
scattering cross-section, which only depends on p and -d for an isotropic interaction potential. Inserting 
Eq. (|A.25p and Eq. (|A.24p in Eq. (jA.23p and performing the angular integral one obtains after lengthy but 
straightforward calculations 



aj(P) = -^^2!^i.EQ (P)^ I'dppV/J (P) f ddsmd{l-coBd)a{d-p) 



(A.26) 
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i=l 



Here, an integration by parts has been exploited. 



d 

dpp'^fip (p) 



dp V - — / ip)a{§;p). 
\ P0 / 



(A.27) 



As a last step, we return to an equation for the original distribution / (P) = ueq (P) X (P)- We recall that 
the equilibrium solution z/eq (P), as defined in Eq. (|4.7p . is a null eigenvector of the Fokker-Planck operator 



^«,<-(P)) + jaj^9(P) 



1 = 1 

such that in particular 

3 

D[7.EQ (P)X(P)]=-^EQ (P)^ 

i=l 

We are thus left with 



(A.28) 
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(A.29) 
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(A.30) 



Switching to dimcnsionlcss variables u = p/pp and using Eq. (j2.6p for (p), one finally has 



3 r 



d M 

(P./(P)) + -^/(P) 



dP, 



[3 dPf 



(A.31) 



with 



16 _ TO 



d'Hw^'e 



di?sini? (1 — cosi?) a upp) , 



(A.32) 



as in Eq. 



Appendix B. List of symbols 



/9 


inverse temperature, /3 = 1/ {kBT) 


Dpp 


diffusion coefficient, see ()5.16p 


Dxx 


diffusion coefficient, see (|5.17p 


V 


friction coefficient, see (]4.33p. (|5.13p 


erf(x) 


error fmiction (198j 


E{Q,P) 


enerffv transfer to test particle, see (|4.1ip 




characteristic function of momentum transfer distribution, see (|7.12p. (]7.15p 


/(P) 


distribution function of the test particle momentum 


/(p/>pO 


elastic scattering amplitude 


/o(p) 


forward scattering amplitude, /o (p) = / (p, p) 


/b(Q) 


scattering amrjlitudc in Born aoDroximation. sec (]2.15|) 




confluent hvpereeometric function |198j 


r(x) 


gamma function |198j 


Ttot 


total scattering rate, see (|7.4p 


r/3 


thermal scattering rate, see (]4.28p 


f(C/) 


dimcnsionlcss loss rate, see (|4.54p 


Ho 


kinetic energy operator of test particle, Hq = P^/ (2Af) 




gas induced energy shift, see (|2.5p. (|6.2p 


K 


scaled momentum transfer, see (]4.52p 




thermal de Broglie wave length of test particle, Ath = ^^271^^ (] / M 


i(p,P;Q) 


non-unitary part of Lindblad operator, see (12.31) 
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Lb (P; Q) non- unitary part of Lindblad operator in Born approximation, see (|2.17p 

C superoperator of the dissipativc part of the quantum linear Boltzmann equation, see (|2.1|) 

fi (p) distribution function of gas momenta 

IJf3 (p) Maxwell-Boltzmann distribution, see (|2.6p ; for n^^'^'^'^^ see (|8.2p . (|8.3p 

m mass of gas particle 

m» reduced mass, ~ mM/ (m + AI) 

M mass of test particle 

(P; Q) classical gain rate, see (|2.11l) 

M^it (P) classical loss rate, see (|2.12p 

superoperator of the quantum linear Boltzmann equation, see ()2.ip 

rigas density of gas particles 

n = ni + in2 index of refraction, see (|6.14p 

i^EQ (P) Maxwell-Boltzmann distribution for the test particle, see (|4.7p 

Pfj most probable momentum in thermal gas distribution, pp — ^^mj (5 

P = (Pi, P2, P3) momentum operator of test particle 

Q momentum transfer to test particle 

set of all momenta perpendicular to Q, i.e., = {p G ^ : P • Q = O} 

p quantum state of motion of test particle 

rel (p, P) relative momentum, see (|2.4p 

(t(P/,P,;) differential cross-section, (T(p/,Pi) = |/(p/,Pi)l^ 

(J differential cross-section (isotropic potential), co^-d ~ pi • Y*} /pf^ Pi — 

(Ttot total elastic scattering cross-section, see (|6.4p . (|7.20p 

(7b (Q) differential cross-section in Born approximation, ub (Q) = |/b (Q)|^ 

S (Q, E) dynamic structure factor, see (|8.12p 

U scaled momentum, see (j4.5ip 

Vfi most probable velocity in thermal gas distribution, vp ~ pf^/m 

V (X — x) interaction potential between test and gas particle 

Kpt optical potential, see ()6.8p 

X = (Xi, X2, X3) position operator of test particle 

X||Q vector component parallel to Q, i.e., X||q = (X • Q) Q/Q^ 

Xj_Q vector component orthogonal to Q, i.e., Xj_q = X — Xyq 
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